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Abstract 

We consider the construction of the deformed Whitham system for the KdV- 
equation in the one-phase case and investigate the conservation of the Hamiltonian 
properties in this situation. It is shown then, that both the Gardner - Zakharov - 
Faddeev and the Magri brackets give the deformed Dubrovin - Novikov brackets (the 
brackets of Dubrovin - Zhang type) for the deformed Whitham system constructed 
by our procedure. The general approach used in the paper gives a scheme for the 
averaging of the Poisson structures in the general situation. 

1 Introduction. 

We are going to consider the conservation of local field-theoretical Hamiltonian structures 
in the method of deformations of the Whitham systems. As it is well known, the Whitham 
method is connected with the slow modulations of parameters of (one-phase or multi- 
phase) periodic or quasi-periodic solutions of PDE's while the Whitham system itself 
rules the behavior of the modulated parameters as of the functions of time and spatial 
variables. The Whitham system is usually written as a system of Hydrodynamic type 



xr T = v»(u)u» (i.i) 

and gives the main term in the connection of the time and spatial derivatives of parameters 
U I/ (X,T). The variables T and X represent usually the "slow" time and spatial variables 
T = et, X = ex connected with the variables t and x through the small parameter 
e. The Whitham system ( 11. ip is then a homogeneous system of Hydrodynamic Type 
connecting the first derivatives of the slow modulated parameters. Different aspects and 
numerous applications of the Whitham method were studied in many different works and 
the Whitham method is considered now as one of the classical methods of investigation 
of non-linear systems. 
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Different properties of the Whitham equations were investigated by many authors (see 
f.i. [I] - 0, © - [13], [21], [29] - [M], [36] - @6], [50] - [52], [M] - [56], [S9] - [M], [66] - [TJ). 

Thus, it pointed out by G. Whitham ([70j EJ E2]) that the Whitham system flU]) has 
a local Lagrangian structure in the case when the initial system has a local Lagrangian 
structure 



on the initial phase space the space {(p(x,t)}. 

The procedure of construction of the Lagrangian formalism for the Whitham system 
(11. ip is given by the averaging of the Lagrangian function £ on the family of m-phase 
solutions of the initial system. Let us note also that in the case of presence of additional 
parameters n l the additional method of the Whitham pseudo-phases should be used. 

The important procedure of the averaging of local field-theoretical Hamiltonian struc- 
tures was suggested by B.A. Dubrovin and S.P. Novikov ([TTJ [T2J EES])- The Dubrovin- 
Novikov procedure gives the local field-theoretical Hamiltonian formalism for the Whitham 
system (II. ip in the case when the initial system has a local Hamiltonian formalism of gen- 
eral type. The Dubrovin-Novikov bracket for the Whitham system has a general form 



and was called the local Poisson bracket of Hydrodynamic type. The theory of the brackets 
(II. 2p is closely related with differential geometry ([HI [121 03]) and is connected with 
different coordinate systems in the (pseudo) Euclidean spaces. Let us say also that during 
the last years the important weakly-nonlocal generalizations of Dubrovin-Novikov brackets 
(Mokhov-Ferapontov bracket and Ferapontov brackets) were introduced and studied ([571 
[23 [2611271 Ml EQl EB [53]). 

During the last years the theory of deformations of systems (II. ip and the Poisson 
brackets (Q was intensively studied ([HIISIIiainitiatlHlI^nill^ 
The e- deformations of systems of Hydrodynamic Type (II .ip and of brackets (11.21) give 
the "dispersive" corrections to (II. ip and (II. 2p and are represented usually as the formal 
series in the powers of e with the higher derivatives of the parameters U. Let us say that 
the theory of the compatible Poisson brackets (11.211 and their deformations demonstrate 
very nontrivial structures and is considered now as one of the general approaches in the 
classification of integrable hierarchies. 

We will consider here the deformations of systems (II. ip and the Poisson brackets (ll.2p 
connected immediately with the Whitham method for the slow-modulated parameters. 
As far as we know, the idea of consideration of dispersive Whitham systems appeared 
first in the paper of M.J. Ablowitz and D.J. Benney (pQ) where the dispersive character 
of the higher corrections in Whitham approach was pointed out. The regular procedure 
of deformation of the Whitham systems was constructed in |55j in connection with the 
theory of deformations of systems of Hydrodynamic Type developed in [T9~l [20] . In [56] a 
special modification of the deformation procedure which gives a regular transition from 
the linear to non-linear systems was also suggested. The procedure used in [56] represents 








{U»(X), U"(Y)} 



g Ufl (XJ) 8\X - Y) + &7(U) U\ S(X - Y) 



(1.2) 
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the method of A.C. Newell giving the connection between the Whitham approach and 
the Nonlinear Shrodinger equation approach to the slow modulations ([58J) which was 
modified to the case of the deformed Whitham systems. 

We will consider in this paper the modification of the deformation procedure considered 
in ([SB]) since the regular properties of the deformation procedure in the case of the 
vanishing amplitude of oscillations seem to be important in many situations. 

The main goal of this paper is to prove the conservation of the local field-theoretical 
Hamiltonian structures in the method of the deformations of the Whitham systems which 
is considered in the example of the one-phase modulated solutions of the KdV equation. 
Namely, we suggest here a scheme of the "averaging" of local field-theoretical Poisson 
brackets giving the deformed Dubrovin - Novikov brackets for the deformed Whitham 
systems (11.11) . The procedure considered here is based on the Dirac procedure of the 
restriction of the Poisson bracket on a sub-manifold which provides the Jacobi identity 
for the "averaged" Poisson structures. 

In Chapters 2 and 3 we describe the scheme of deformation of the Whitham system for 
the KdV equation giving the dispersive corrections to the standard system of Whitham 
in this situation. In Chapter 4 we consider immediately the averaging of two local Hamil- 
tonian structures for KdV and prove the existence of two deformed brackets fll.2|) for 
the deformed Whitham system. Finally, in Chapter 5, a scheme of the averaging of the 
local Lagrangian structures in the method of deformations of the Whitham system is also 
considered. The construction used here have in fact a general character and can be used 
in analogous form for different systems of PDE's. 

2 The Whitham method and the deformation scheme. 

As is well known the Whitham method ([701 ED E2]) is connected with the slow modula- 
tions of periodic or quasiperiodic m-phase solutions of nonlinear systems 

F\cp,<p t ,<p x ,...) = , i = l,...,n , <p = {<p 1 ,...,<p n ) (2.3) 
which are represented usually in the form 

^(x,t) = $*(k(U)x + w(U)* + 0o, U) (2.4) 

In these notations the functions k(U) and w(U) play the role of the "wave numbers" 
and "frequencies" of m-phase solutions and 6q are the initial phase shifts. The parameters 
of the solutions U = (U 1 , . . . , U N ) can be chosen in arbitrary way, however, we assume 
that they do not change under arbitrary shifts of the initial phases 6 Q of solutions. 

The functions $ l (0) satisfy the system 

F i (<S>,cu a <S>8°,kP<S> e P,...) = , i = l,...,n (2.5) 

and we choose for every U some function 3>(#, U) as having "zero initial phase shifts" . The 
full set of m-phase solutions of (12. 3p can then be represented in form (12.41) . For m-phase 
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solutions of Q we have then k(U) = (V(U), . . . , k m (U)), w(U) = (^(U), . . . , w m (U)), 
#o = • • • , # m ), where U = (U 1 , . . . , U N ) are parameters of the solution. We require 
also that all the functions $*(#, U) are 27r-periodic with respect to every 9 a , a = 1, . . . , m. 

Let us denote by A the family of the functions 3>(#, U) which depend on the parameters 
U in a smooth way and satisfy system (12. 5p for all U. We will assume also that A is the 
maximal family having these properties. 

In Whitham approach the parameters U become slow functions of x and t: U = 
U(X, T), where X = ex, T = et (e -> 0). 

The functions U(X, T) should satisfy in this case some system of differential equations 
(Whitham system) which makes possible the construction of the corresponding asymptotic 
solution. More precisely (see [50]), we try to find the asymptotic solutions 

<p\0,X,T) = Y,&w(^^+0,X,t) e k (2.6) 

fc>0 ^ ' ' 

(where all VE^) are 27r-periodic in 0) which satisfy the system (12.31) . i.e. 

F l (ip,eip T ,ecp x ,...) = , i = l,...,n 

The function S(X,T) = (S X (X, T), . . . , S m (X, T)) is called a "modulated phase" of 
solution (1231) . 

It is easy to see that the function ^/^(0, X, T) should belong to the family of m-phase 
solutions of (12. 3p at every X and T. We have then 

* (o) (0,X,T) = $(0 + o (X,T),U(X,T)) (2.7) 

and 



S%{X,T) = cu a (XJ) , S%(X,T) = k a {U) 

as follows from the substitution of 1 12. 6ft into system ( 12. 3ft . 

The functions ^ r ( / t)(0, X, T) are defined from the linear systems 

V j[XJM {X,T)¥ {k) {e,X,T) = f {k) {0,X,T) (2.8) 

where L % jm g i (X, T) is a linear operator given by the linearization of system (12. 5p on solu- 
tion (12. 7p . The resolvability conditions of systems (12.81) can be written as the orthogonality 
conditions of the functions f(jt)(0, X, T) to all the "left eigen vectors" (the eigen vectors of 
adjoint operator) of the operator L 1 ^ g AX,T) corresponding to zero eigen- values. The 
resolvability conditions of (12. 8p for k — 1 

V mM {X,T)¥ {1) {e,X,T) = f {l) {0,X,T) (2.9) 

together with the relations kj. = give the Whitham system for m-phase solutions of 
( 12. 3p which plays the central role in the slow modulations approach. 
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Let us say that the resolvability conditions of (I2.8P can in fact be rather complicated 
in a general multi-phase case. Indeed, we need to investigate the eigen-spaces of the 
operators L[u,e ] and LiLj g , on the space of 27r-periodic functions which can be rather 

non-trivial in the multi-phase situation. Thus even the dimensions of kernels of -£/[u,0 o ] 
and LLj can depend in non-smooth way on the values of U so we can have a rather 
complicated picture on the U-space ( [HI LTl [H] ) - 

These difficulties do not usually appear in the one-phase situation (m = 1) where the 
behavior of eigen-values of £[u,0 o ] an d ^lu e ] * s usuan y rather regular. 

In this chapter we are going to consider a scheme of deformation of the Whitham sys- 
tem giving "dispersive" corrections to the system of Hydrodynamic Type which describe 
in fact the higher corrections to the corresponding asymptotic solutions. We are going 
to use here the one-phase modulated solutions of the KdV equation as a basic example 
throughout the paper, so let us consider now the KdV equation 

ipt + <p<p x + ip xxx = (2.10) 
It has a family of exact solutions of the form 

<p(x,t) = §(kx + ut;U\U 2 ,U 3 ) (2.11) 

where the functions $(#; U 1 , U 2 , U 3 ) depending on three real parameters U 1 , U 2 , U 3 is 2tt- 
periodic in 9, the wave number k and the frequency u are uniquely determined by these 
parameters, 

k = k(U\U 2 ,U 3 ) , u = lu(U\U 2 ,U 3 ) 

As we pointed out already the Whitham modulation theory gives a prescription for finding 
approximate solutions to KdV in the form 

$ ^ S ^ T) ; U l (X, T) , U 2 (X, T) , U 3 (X, T)^j (2.12) 

where e is a small parameter, 

X = ex , T = et 

are slow variables, the dependence of the parameters U u = U U (X,T) is determined from 
certain hyperbolic system of the first order quasilinear equations of the form 

U» = V;(U\U 2 ,U 3 ) U£ , u,fi = 1,...,3 (2.13) 

The phase function S(X, T) is determined by quadratures 

d x S(X,T) = k(U\X,T),U 2 (X,T),U 3 (X,T)) 

(2.14) 

d T S(X, T) = u(U\X, T), U 2 (X, T), U 3 (X, T)) 

We now want to describe the higher corrections to the approximate solutions (12.121) . We 
will call them deformed Whitham equations. 
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The deformed Whitham equations will arise in the description of solutions to (I2.10p 
in the form 



p = <S>(S(X,T) + 8;U(X,T)) + J2$(i)(S(X,T) + 6;X,T) (2.15) 

i>i 

where the functions 

a>®(6)X,T) = $^;U,u x ,u xx ,...,u«) 

27r-periodic in 8 are graded homogeneous differential polynomials in Uj, XJxx etc. with 
coefficients being smooth functions of U = (U 1 , U 2 , U 3 ). The gradation is defined by the 
rule 

degd^U = m, m = 1,2,.... 

As usual the degree of the product of homogeneous differential polynomials is equal to 
the sum of their degrees. We use here the notations X and T just to emphasize that the 
functions \J(X,T) are "slow" functions of spatial and time variables. At the moment we 
do not write the small parameter e explicitly; it will be reintroduced later. 

It will be convenient to choose a particular system of coordinates U ,U 2 ,U 3 in the 
space of traveling wave solutions U). We denote them 

U = (k, co, n) 

where k and u are the wave number and the frequency and n is the mean value of $. The 
ODE for the function <3> 

cu$ e + k^ e + k 3 ®eee = (2.16) 
can be integrated by quadratures 

A? f d§ 

q = 




2 J as v/-£;$ 3 /6 - w$ 2 /2 + A$ + B 



where 03 is the third zero of the cubic polynomial — /c$ 3 /6 — uj§ 2 /2 + + B according to 
the normalization shown at Fig. [TJ The dependence on the parameters of the coefficients 
of the polynomial A = A(k,u,n) B = B(k,u,n) is determined from the equations 



k 3 f rf$ 



2 J ^/-k<5> 3 /6 - u<5> 2 /2 + A& + B 
k 3 f $6?$ 




2tt 



2nn. 



2 / y/-k$ 3 /6 - u$ 2 /2 + A& + B 

We also fix the initial phase shift of the functions $(#, k,u,n) in such a way that every 
$(#, k, u, n) has a local maximum at the point 8 = (see Fig. [T|). 

It is well-known that the function — $(A;x + ut, U)/6 represents the one-gap potential 
for the Shrodinger operator 
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Figure 1: The function $(#, k,u) having zero initial phase shift. 



L = 



d (p 
dx 2 6 

while the KdV equation can be written in the Lax representation 



(2.17) 



where 



dL 
~dt 

dx 3 



A, L 



^Tx + Tx* 



Let us say that the integrability of the KdV equation will be convenient in some aspects 
of our considerations. However, the general questions considered here are not connected 
with the integrability and are applicable for the non-integrable examples as well. 

It is well-known also that the solution $(kx + cut, U) can be represented in the form 



§(kx + ut, U) 



2a 



dn 2 



+ 7 



V 



2a 
3? 



(2- S 2 ) + 7 



where s is the modulus of the Jacobi elliptic function dn(w,s), < s < 1. The value 
2a plays the role of the amplitude of oscillations for the periodic solution and the values 
(k, cj, n) can be expressed in terms of the parameters (a, s, 7) in the following way 
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7T / a \V2 4-7T 2 . / a \ 3 / 2 77T / a \ V2 

n = 7 + 



K(s) \Qs 2 J ' K(s) v y V6s 2 / K(s)\Qs 2 

2aE(s) 



s 2 K(s) 

where K(s) and E(s) are the elliptic integrals of the first and the second kind respectively. 
We can write also 

*(M,*,7) = |dn 2 (^, S ) + 7 

for our normalization of the functions U). Let us note also that the parameters 
(a, s, 7) are connected with the energy band edges (r 1; r 2 , r 3 ) of the operator (I2.17P by the 
formulas 

r-2 - n 2 

r 2 - ri = a , — s , n + r 2 - r 3 = 7 

r 3 - r x 

(r 3 > r 2 > ri). 

The total function $( to *)(0, X, T): 

$( to ')(^,X,T) = ^T$ (0 (fl,X,T) = <l>(e-S(X,T),X,T) 

l>0 

satisfies the equation 

+ 35 x $^ + 3S XX $^ } + S xxx ^ ot) + (2.18) 

This yields linear equations for the functions u) for I > 1. In particular the function 
X, T) satisfies the equation 

w$(i )fl + fc$ ( i)$e + + A; 3 $ (1)ee9 = f {1) (9,X,T) (2.19) 

where 

f {1) (9,X,T) = " - 3A; 2 $, ex - 3A;£; x $, e (2.20) 

and the notation M means that we consider just the part of <&t having degree 1 according 
to our definition. 

Denote L\x,t\ the linear operator 
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We can rewrite (I2.20p in the form 

L[X,T] $(1) = /(l) 

In the same way we have the analogous systems for the functions $m(0, X, T) having the 
form 

L [x ,T]^(i) = u$ ( i )e + fc$(i)$e + fc$(j)*$ + fc 3 $(z)ew = f (l) (9,X,T) (2.22) 

where fm(9, X,T) are the discrepancies having degree /. 

The functions k(X,T) = Sx, uj(X,T) = St and n(X, T) must satisfy the "deformed 
Whitham system" 



^ <7(q (fc, u, n, k x , Ux, %,...) (2.23) 

Z>1 

n T = y^77 ( ;) (k,u,n,k x ,wx,n x , ■ ■ ■) 

2>1 

where all cm, 77/n are graded homogeneous differential polynomials in (k, co, n, k x , uj x , n x , ■ ■ 
of the degree I. 

It is easy to see that relations f!2.23j) give in fact a possibility to represent in the form 
of homogeneous differential polynomials any expression like krx...x, ^tx...x, nrx...jf, and 
even hr...TX...Xi w T...rx...x, n T...rx...x iterating the subsequent substitution of the series 
(I2.23p . (The last property will not be necessary for the KdV equation). 

According to (I2.23P all the time derivatives like $t, $(z)t can also be represented as 
the sum of homogeneous components 

- $ (Z)T + $ (Z)T + $ (Z)T + • • • 

where the functions $[n T are differential polynomials of (k, to, n, kx, ^x, nx, ■ ■ ■ ) of the 
degree s. 

We impose the following orthogonality conditions on the discrepancies f(i)(0, X,T) 
2t , dO f 2n dO . 

f<» 2i = ■ | * / «" S = (2 ' 24) 

and also the "normalization" conditions 



* (0 ^ = , ^ * w - = (2.25) 

for the functions T) defined from (I2.22p modulo the linear combinations 

a(X,T)$ e + b(X,T)$ n . 
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For determination of cm, tjq) we use the system f!2.23j) to remove all time derivatives 
of (k, u, n) after the substitution of (12.151) into (12.101) in order to represent (I2.18P in the 
graded form. 

The functions a^, r]m arising in (I2.23P are found from the compatibility conditions 
of systems ( I2.22p in the Z-th order. It can be shown that conditions (I2.24l) - (l2.25p define 
uniquely all the expressions am, r/^ and the corrections Qm, I > 1. 

So, our prescription for deriving the system f !2.23|) is based on the following three 
conditions: 

I) All the functions k,u,n) are chosen in the way shown at Fig. 1 ; 

II) The modulated phase S(X,T) is connected with the parameters (k,u,n) by the 
relations 

S T (X,T) = u(X,T) , S X (X,T) = k(X,T). 

III) All the higher corrections §(i)(9, X, T), / > 1 satisfy the normalization conditions 
(I2T25D . 

According to the statements above system (I2.23|) is uniquely defined by the conditions 
(I)-(III). 

Let us now say some words about solutions of the system (I2.22p . We have 

+ + 
Put $(q = to arrive at 

So 

k 3 a { i )e ($e) 2 = 

and 

(2.27) 

Formula (I2.27P has a local character and we have to investigate solution (I2.27P on the 
whole axis — oo < 9 < +oo. (The expression l/($g) 2 has singularities at the points 
9 n = 7m, n G Z). Two important cases can be pointed out in our situation: 



(1)60 



J f {l) {9',X,T)d9' + 6 (2.26) 



2k 3 a {l )8<S> ee = I f(i)d9' + 6- 



®o> I f(i) d9" + 6 



d& + ^ 



e i re 



l — J $ fl „ J f {l) d9'"d9"d9' 
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I) The function f(i)(9) is even, f(i){-9) = f(i){9); 

II) The function f (l) (9) is odd, f {l) (-9) = -f Q) (9). 

Proofs of the following two propositions are straightforward. 



Proposition 2.1. 

For an even smooth periodic /(«(#) the corresponding solution $>m(9) of ( QL 22\) satis- 
fying conditions 112.25]) is an odd smooth periodic function. 

Proposition 2.2. 

For an odd smooth periodic fm (9) the corresponding solution <£>(/) (9) of Ii2.22\) satisfying 
conditions \2.25S) is an even smooth periodic function. 

In particular the function f^(9,X,T) is given by the following expression 

— [ < £(o)t]^ — ®(o)®(o)x — 3Sx$(o)eex — 3S x S X x®(o)ee- 

Recall that the mark M means that we collect the terms of the degree 1 . 
Orthogonality conditions (12. 24ft can be written as 



»2tT 



2tt 



^ d9 

*(0)n— V(l) + 



2k 



d9 



% dB 



and 



I 



2w ~ d9 f 27T x d9 f 2n x d9 f 2n ,. dO 

^(O^Wco— (7(1) + J ^(0)$(0)n— V(l) + J ®(0)<5>(0)k— UJ X = J f {l) — 



where 

/(i) = — ®(o)®(o)x — 3S x §( )ggx — 3S x S X x^{o)ee 

The equations written above determine the functions o"(i)(/c, u, n, k x , u x , n x ) and 
77(1) (/c, u, n, k x , u x , n x ). In this way we arrive at the standard Whitham system (12.131) as 
the zero order approximation of (I2.23p . 

More generally, according to our approach the functions f^(9,X,T) will be always 
represented in the form 

f(l) = ~ [$(0)T] + /(l) = - $(0) W <T(i) + $(0)n V(l) + /(I) 

where does not contain the terms <j(;) , rj^ . The corresponding orthogonality conditions 
(12.241) recursively determine all the terms cr(j), t/(z). 

It is easy to see that the function /(i) is even: f(i)(—9) = f(i)(9). We obtain therefore 
that the function $(i)(0) is odd $(i)(-0) = -$ ( i)(6>). 
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Furthermore, a direct substitution gives 



/(2) 



^(0)^(2) + ®(.Q)nV(2) + f(2) 



where fi 2 \ (9) is odd. 

Using equations (I2.24p for I = 2 we get immediately ct( 2 ) = 0, 77(2) = for the next 
terms in f 1 2 . 2 3 1) . The total function f^) (9, X, T) becomes then an odd function in 9. Hence 
the second correction $(2)(#) is even. By simple induction we obtain the following Lemma: 

Lemma 2.1. 

For the choice of the functions $(#; k,u,n) corresponding to Fig. [1\ the following 
statements are true: 

1) All the even terms a^i) (k,u,n, . . .) , rj^i) (k,u,n, . . .) in the deformation of Whitham 
system \2.23\) are identically zero: a^i) = 0, r/^i) = 0; 

2) All odd corrections ^(2i+i)(9,X,T), I > in I12.15\) are odd functions in 6; 

3) All even corrections §(2i){9, X,T) , I > 1 in A2.15\) are even functions in 9. 

3 Deformation scheme for the case of small ampli- 
tude oscillations. 

The above procedure of deformation has one weak point. Namely, in the procedure 
described the higher corrections ^>^(9,X, T) as well as the higher deformation terms in 
system (I2.23P are singular in the limit of small amplitude oscillations of ip(x,t). The 
reason for such a singular behavior can be explained in the following way. 
Let us rewrite system (12.221) in form (12.261) i.e. 



+ fc$ (0 $ + k 3 <5> {l)ee = g { i){9,X,T) 



where the right-hand part g^ given by the expression 





is periodic in 9 due to the conditions (12.241) . 
We can rewrite this system in the form 



Q{X,T]®(1) = 9(1) 



(3.1) 



where 



Q [x ,t\ = u(k,A,n) + k<$> + k 3 — 



(3.2) 



is a self-adjoint operator on the space of 27r-periodic functions. 
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Operator (13.21) has just one eigen-vector $0 with zero eigenvalue on the space of 2ir- 
periodic functions. The constants 8(i)(X,T) are uniquely determined by the second con- 
dition (I2.25P for the solutions of (13. ip . It is not difficult to get analytic expressions 
for 6(i)(X, T). It is also easy to see that Sq)(X, T) = for I = 2s + 1, s > 0. 

Provided that conditions (I2.24p are satisfied we can write the solution of ( 13. ip in the 
form 

(3.3) 



(0 



Ex-^' X > T ) fag®) 



where £j(9,X,T) are the normalized eigen- vectors of Q\x,T] corresponding to non-zero 
eigenvalues Xj. 

Let us consider now operator (13. 2p for the case of small amplitude oscillations: 



$(0,X,T) = n{X,T) + a {X,T) cos 9 



a — > 



The parameter ao(X, T) is the amplitude of the first Fourier harmonic of $(#, X, T) which 
is similar to the parameter A = $ max — $ m j n in the limit A — > 0. 

Operator (13.21) has always the eigen-vector £(9,X,T) = $g(9, X,T) corresponding 
to zero eigenvalue, which corresponds to the function — sin 9 in the limit A — > 0. The 
dispersion relation u = u(k,A,n) becomes the dispersion relation of the linear system 



00 



-nk + k 



for A = 0. 

However, the function cos 9 gives also an eigen-vector of linear operator (A = 0) (13.21) 
corresponding to zero eigenvalue. As a result there exists an eigen-vector £i(9,X,T) of 
the operator Q\x,t] corresponding to "small" eigenvalue Ai — > (for A—*0). 

The values and Ai can be also expressed in terms of the elliptic functions in 

our case. Indeed, if we compare the operator —Q[x,T\/k with the Shrodinger operator 
(I2.17P we can easily see that the operator —Q[x,T]/k is represented by the Shrodinger 
operator with a one- zone potential multiplied by 6. It is well known (see f.e. [3 0]) that 
the operator —Q[x,T]/k represents a 3 - energy gaps Shrodinger operator with an elliptic 
potential in this case. 

The spectrum of the operator —Q[x,T]/k is shown at Fig. [2] and we are interested here 
in the particular functions £i(#) and Ai. 

It's not difficult to see then that the eigen-function <&e(9) and £i(#) correspond to 
the gap edges E% = and E4 = —\\/k. Easy to see also that in the limit of the 
small amplitude of oscillations we have \E± — E%\ — > in the full accordance with the 
perturbations theory. The expressions for the functions and Ai can be written in 

the form (see El): 



£i(0, a, s, 7) 



dn 



K(s) 



71 



2s 2 



vT 



As 4 



5s sn 



K(s) 



71 



e,s 
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Figure 2: The spectrum of the operator —Q[x,T\/k. The intervals [E ,Ex] , [E 2 ,E^\, 
[E4, E 5 ] and [Ex, E 2 ], [E 3 , E4], [E 5 , E 6 ] represent the energy bands and the energy gaps of 
a finite size respectively. 



3/^2 

7T 



Ai(a,s,7) = - K 2 (s) (2VI -s 2 + 4s 4 - 2 + s 2 ) k 

= _ n ^2Vl-s 2 + 4s 4 - 2 + /) (^) 3/ y #00 

in our notations. 

By direct substitution it is not difficult to also get the following relations for the values 
of cj(k, A, n), $, Q, £1 and Ai: 

tu = -kn + k 3 -^ j - + O(a 4 ) (3.4) 

a 2 

$(0,k,A,n) = n + a cosfl + — ^ cos 20 + 0(ajj) (3.5) 

^2 ^2 ^2 

Q[M,»] = ^ 3 - ^ + cos^ + ^ cos 20 + k 3 — + 0{al) (3.6) 
U0,k,A,n) = cos6 - £jL + £jL cos2 6 + O(a 2 ) (3.7) 

Ar = -g| + O(aJ) (3.8) 

We can see that the solutions ( 13.31) become singular in the limit of the small amplitude 
oscillations A — > if we do not put additional requirement 

(£1, <?(/)) = 

for all gyy 

The idea of correction of the Whitham approach for the almost linear case using a 
correction of a dispersion relation was invented by A.C. Newell (see [5B], Chapter 2). 
In [5B] the method of A.C. Newell was generalized to the case of the deformed Whitham 
systems which gives the regular deformation procedure in the limit of the small amplitude 
of oscillations. 

To improve the deformation procedure described above we will use the deformation 
scheme suggested in for the case of almost linear systems. 
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Namely, the orthogonality of g(i)(0, X, T) to X, T) can be provided in the following 
way: 

First of all we choose the parameters (k,A,n) instead of (k,u,n) as the regular pa- 
rameters everywhere (including the region A — > 0). Now the main approximation in the 
asymptotic solution ( I2.15P will be again given by the function $>(S(X,T) + 9, k,A,n) 
such that Sx{X,T) = k(X,T). So we have again the same approximation with the same 
relation between S and k as previously at every T. However, we make now also the "de- 
formation" of time evolution of phase S(X,T) such that St{X,T) ^ u(k,A,n) anymore. 
Instead, we put now the deformed relation 

S T = u(k,A,n) + ^u^k, A,n,k x , A x ,n x , ■ ■ ■) (3.9) 
i>\ 

connecting the time derivative St and the parameters (k, A, n) of the main approximation. 
Here again all the functions ujn\{k, A, n, kx, Ax, nx, ■ ■ ■ ) are differential polynomials in 
(kx, Ax, n>x, ■ ■ ■ ) of the degree / with coefficients smooth in (k,A,n) according to the 
same gradation rule, i.e. 

all the functions f(k,A,n) have degree 0; 

the derivatives kix, A [X , rii X have degree Z; 

the degree of the product of homogeneous differential polynomials is equal to the sum 
of their degrees. 

As we have already said the parameter A = § max — Q m i n plays here the role of the 
amplitude of oscillations and we have A(X,T) ~ a (X,T) for the small A. 
We write now the deformed Whitham system in the form 

k T = I w(k,A,n) + y^U(i)(k,A,n,kx,Ax,n x , ■ ■ ■) 
V i>i 

A T = / u at® (k, A, n, k x , A x ,n x ,---) (3.10) 

i>i 

n T = ^ n(f) {k, A, n,k x ,A x ,n x ,---) 
i>i 

which gives a full deformation of the Whitham system having a regular behavior in the 
case of small amplitudes. 

The functions am, rjm are defined as previously from the orthogonality conditions 
of the functions f(i)(&, X,T) to the "left" eigen- vectors $(#) and 1 of the operator L 
corresponding to zero eigenvalues. The functions oum in (13 .9p are defined now from the 
orthogonality of the functions gm(9, X,T) to the eigen-vector £i(9,X, T) of the operator 
Q[x,t] corresponding to the "small" eigenvalue Xi(k,A,n). 

So now we have the condition 

/27T jn 
b(9,X,T)g {l) (9,X,T)— = (3.11) 
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in addition to conditions (I2.24p . The functions \\(k, A, n), k,A, n) are denned by 
continuity on the whole family of one-phase solutions so we can define the system (I3.10P 
on the full space of parameters. 

For our choice of the functions k,A, n) it is easy to prove that the function 
£i(6>, k, A, n) is even in 9. 

For the solutions Q^(9,X, T) we will have automatically 

£i(9, X, T) &q)(9, X, T) — = (3.12) 

in addition to normalization conditions (I2.25p . 

In the same way as previously the following lemma can be proved for systems ( 13 .9 j) - 
(13.1 Op and the asymptotic expansion 

</>{6,X,T) = ®(S(X,T) + 0,k,A,n) + y $2$ ( n(S(X,T) + 6,X,T) (3.13) 

z>i 



Lemma 3.1. For the "unified" choice of the functions <&(9,k,A,n) corresponding to 
Fig. U\ the following statements are true: 

1) All even terms a^i) (k,A,n,...), rjpi) (k, A,n, . . .) in the deformation of the Whitham 
system A3. 10}) are identically zero: a/^i) = 0, r/^i) = 0; 

2) All odd terms u)(2i+i)[k, A, n, . . . ), I > in the deformation A3. 9}) of the dispersion 
relation are identically zero: W(2«+i) = 0; 

3) All odd corrections <&(2i+i){9,X,T), I > in A2.15}) are odd in 9; 

4) All even corrections ^^i){9, X,T) , I > 1 in A2.15}) are even in 9. 

So we can rewrite the relation (13. 9p and the system (13.101) in the form 

S T = u(k,A,n) + ^2uj(2i){k,A,n,kx,Ax,n x ,...) 
i>i 

k T = iu(k,A,n) + y^U(2i)(k, A,n, k x , A x ,n x , ■ ■ ■) 

V i>i 

(3.14) 

A T = ^ of( 2t+ i)(fc, A, n, k x , A x , n x , . . . ) 

l>0 

n T = ^ V(2i+i) (k, A,n,k x ,A x ,n x ,...). 

l>0 

We can see that for our choice of the functions $(#; k, A, n) the full deformation (13.141) 
of the Whitham system includes only odd degrees of the expansion in higher derivatives 
which emphasizes the dispersive character of the deformation. 



V 
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To calculate the terms up), «(3), V(3) ^ us write down the expressions for the discrep- 
ancies /(i), /( 2 ), /( 3 ) in the form 

-/(i) = $? ] + ^x + 3S x S xx $ee + 3S x $ eex + S [ ^ e 

-f(2) = $t + $[ 2 i )T + 3S xx <S> ex + 35 X $ 9XX + + SV$(i)$(i)0 + 

+$$ ( i )x + $ ( i)$x + 3S x S xx $ (1)ee + 3^$ ( i )eex + SP$ e 

-/(3) = $t + $[i) T + $( 3 2 ] )r + ®xxx + 3S xx <5> (1)ex + 3S x $ (1)exx + S xxx $ {1)e + $ ( i)$ ( 
+$$ (2)x + $ (2) $x + 3S x S xx $ (2)ee + 3S x <f> {2) e0x + 5 x $ (1) $ (2)e + S x $ (2) $ ( i )9 + SP$ 9 . 

Let us remind again that we do not prescribe any certain degree to the operator d/dT, 
so we have 

S T = uj(k,A,n) + u^k, A,n, k x , A x ,n x ) + ■■■ 
<£ r = $? ] +4 2] +4 3l + ... 

$ (1)T = + $Si)T + *(i ] )T + - 

cf, _ r&[3] i cBW , *[5] , 
*(2)T - $ (2)T + $ (2)T + * (2)T + • • • 

where all ^j-^ are differential polynomials in (k x , A x , n x , . . . ) of the degree s with smooth 
coefficients depending on (k,A, n). 

It is easy to see that the only odd in 6 term in is — S^Qg. So from the orthogonality 

of g(i) to £i(6,X,T) we get immediately W(i) = Sjr = in accordance to Lemma 2. In 
the same way we put also oj^ = 0, 5^ = for all / = 2s + 1, s > 0, such that only u^ s ), 
dps) should be computed for s > 1. 

We will not need to calculate completely the system (I3.14p in this paper, however, let 
us briefly describe here the scheme for the determination of the functions 77(1), «( 3 ), 
77(3), up). We have 



$J ] = ^ a Af + $ n r^ ] + $ fc 4 1] = $ a« ( i) + $„r/(i) + (u(k, a, n)) 



x 



so the orthogonality of f^(9,X,T) to the functions $(6*,X, T) and 1 gives the usual 
expression for rj^ given by the standard system of Whitham. The only even term 

f2l 

in fp) is — $ T . So we get immediately «( 2 ) = 0, ?7( 2 ) = from the orthogonality of /( 2 ) to 
$(#,X, T) and 1. The term — ®f\ T is given by 

w a(i)( ) + "wr y( ] ) 
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and it is a known function. From the orthogonality of g@) to £1(8, X, T) we get a relation 
for u)(2)(k,A,n, . . . ): 



where 



9(2) 



c2ix in 

^ I 6(0) (HO) - *(o)) - 



$!J T + 3/cx$e'x + 3/e$0/xx + + + $$ ( i)x + 



2tt 



6(0)<?(2)W 



2^ 



+$ ( i)$x + 3fcfc^$ ( i )0 / fl / + 3fc 2 <^ (1)f? ^^] dd' + 5 (2 )(X,T) 

It convenient to determine the values of u^) and 5^) simultaneously from the orthog- 
onality of g(2){0) to both the vectors and 



U(0) 

and add the relation 

W(2) 



^2s £2 



J: 



2.s 



27T Jfl 

u(o) (m - $( )) - 

Zix 



2tt 



2 \ (0) de 



df) 

U(e) g' m {6) - 



which gives a non-degenerate linear system on the values Uf 2 ), 5(2)0 

[31 

Repeating all the arguments we get W( 3 ) = = from the orthogonality of g^ to 
cfl v ^ w e have also $[f )T = view A^ 1 = 0, n§ = 0, A^ 1 = 0- The function -$[ 3 2 ] )r 



is given by 



SA(Y) 



a ( i)(y) + 



<Jn(y) 



<jjfe(y) 



u Y (Y) dY 



and is a known function again. 
We have then 

$T = $A«(3) + $n?7(3) + ( w (2)) x 

so we get the functions a(3), 77(3) from the orthogonality of /(3) to $>(9,X, T) and 1. 

It is easy to see also that the procedure can be extended to any order / such that all 
the terms W(2j), cxpi+i), V(2i+i) wm be uniquely determined. 



The system f)3.14p determines the evolution of the parameters (k, A, n) of the zero 
approximation of ( 12.151) such that the following conditions are satisfied: 

I' ) All the functions k, A, n) are chosen in the way shown at Fig. 1. 

II' ) The modulated phase S(X, T) and the parameters (k, A, n) of the zero approxi- 
mation are connected by the relation 



1 Both the functions £,tot{&) have in fact explicit expressions in terms of elliptic functions. 
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S X {X,T) = k{X,T). 

Ill' ) The higher corrections $(/)(#, X, T), I > 1 satisfy normalization conditions (I2.25P 
and f l3TT2|) . 

We would like to introduce a small parameter e according to our gradation rule for 
more convenient notations. System (I3.14p will be rewritten in the form 

S T = ui(k,A,n) + ^2e 2l u { 2i)(k,A,n,kx,A x ,nx,---) (3.15) 



i>i 



I u;(k, A,n) + ^e 2l i 
\ i>i 



k T = \ u(k,A,n) + \ e u {2 i)(k,A,n,k x ,A x ,n x ,.. 



x 



A T = ^2e 2l a ( 2i +1) (k,A,n,k x ,A x ,n x ,...) (3.16) 
z>o 

n T = ^2e 2l r](2i+i)(k,A,n,k x ,A x ,n x ,...). 

l>0 

The asymptotic expansion (I3.13P will also be rewritten in the form 

X,T) = ®(?^ + 9,k,A,r?) +E e ^(o(^T^ + ^^ T ) ( 3 - 17 ) 



Z>1 

according to the gradation rules for the functions S(X, T) and $(;)(#, X, T). 

In these new notations we can actually see that system f l3.22l) -( l3T23l) describes the 
asymptotic solutions of the equation 

V?t + Wx + £ 2 <fxxx = (3.18) 

where the small dispersion e 2 arises after the rescaling T — > eT, X — > eX. So our further 
considerations will be applied to the KdV equation in the small-dispersion form (I3.18p . 

Let us emphasize, however, that (13.171) is not an e-expansion of the asymptotic solution 
of (I3.18P since all of the functions k(X, T, e), A(X, T, e), n(X,T,e) are solutions of e- 
dependent system (I3.16p . such that expansion (I3.17P can contain more complicated e- 
dependence. According to our rules we should not separate the different orders in e of 
the functions k(X, T, e), A(X, T, e), n(X, T, e) and just use the gradation rules formulated 
above for the e-dependent functions. □ 

The solutions of system (13 . 1 6[) can be considered in different ways. Thus, it is easy to 
define the formal graded form of solutions of ( 13 .161) 



2 As it was pointed out in [55] the series (|3.17|) (or (|3.13p ) corresponds to the expansion with respect 
to X-derivatives of the "renormalized" e-dependent parameters of the main approximation &(8,k,A, n) 
which gives these specific rules of constructing of the series (|3.17[) . 



19 



k(X,T) = k(X,0) + ^T i e / - 1 K (0 (A;(X,0),A(X,0),n(X,0),...) 

i>i 

A(X,T) = A(X,0) + Y,T l e l - 1 A {l) (k(X,0),A(X,0),n(X,0),...) 

l>i 

n{X,T) = n(X,0) + J2T l e l ~ 1 N {l) (k(X,0),A(X,0),n(X,0),...) 



i>i 

for < T < 5, where all K^, A^, are local functionals of k(X, 0), A(X, 0), n(X, 0) 
and their derivatives having the corresponding degree. 

However, a more complicated treatment of the solutions of ( I3.16P connected with their 
global behavior based on the so-called quasitriviality transformations (see [19j 120]) of 
parameters (k, A, n) is also possible and seems to be very important in the theory of the 
deformed Whitham systems. 

Let us recall also that the KdV equation ( 12.101) has an infinite series of conservation 
laws which can be written in the form 

±P»(<p, %,...) = -^-(yfo p.,...) (3.19) 

(Jj V KAjJii 

with some local functionals V u ((p, <p x , . . . ), Q u (f, f x , ■ ■ ■)■ For equation (" 13 . 1 8[) the corre- 
sponding relations can be written respectively 

-^V v { V ,e Vx ,...) = ^Q%,e^,..) , v = 0,1,2,... (3.20) 
According to standard numeration we put 

V° = <p , Q° = -\^-e\ xx 
for v = and we have the conservation of the Casimir function 

/+oo 
ipdX 
-oo 

for the Gardner - Zakharov - Faddeev bracket in this case. 
For v = 1 it is put traditionally 

T*>1 ^ 2 -^3 2 1 2 2 

V = ' Q = ~3 V ~ 6 Wxx 2 6 Vx 

which corresponds to the conservation of the momentum functional for the same bracket. 
For v = 2 we put 

V 2 = \ ^ - g e Vx > Q 2 = - g ¥> 4 - 2 e 2 <P 2( Pxx + (?<P<Px + t A <Px<Pxxx - ^ eVL 
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which gives the conservation of energy in the Gardner - Zakharov - Faddeev Poisson 
structure. 

The higher conservation laws are connected with the integrable nature of the KdV 
equation and arise from the method of the inverse scattering problem. 

It's not difficult to see that the conservation laws of the KdV equation (I3.18P give 
conservation laws for system (I3.16P after the "averaging" of the corresponding densities 
V v , Q u on the asymptotic family (I3.17p . Indeed, after the substitution of (13. 17j) into 
(I3.20p and integration w.r.t. 9 we get the relations 

±<T) = ^<C> (3.2D 

where the quantities < V v >, < Q u > are given by the substitution of solutions (13. 17ft in 
the expressions for V v and Q u and integration w.r.t. 6 over the period. It's not difficult to 
see also that the values < V v >, < Q v > are expressed in this case as the local functionals 
of the parameters (k, A, n) and their X-derivatives 



(V) = {V v ){k,A,n,k x ,A x ,n x ,...) , (Q v ) = {Q v ) {k, A,n,k x ,A x ,n x ,...) 

which are polynomial in the derivatives of (k, A, n) and can be written in the graded form 
we introduced above. Relations (I3.2ip give then an infinite set of conservation laws for 
system (I3.16P written in the same graded form. 

We can see also that the values of independent integrals < V v > can be also chosen 
as the parameters of the solutions (13.1 7ft such that the values (k, A, n) will be expressed 
in the form of graded expansions with respect to the X-derivatives of < V v > given by 
the "inversion" of the corresponding expansions for < V v >. System (I3.16P written in the 
corresponding parameters (say < V° >, < V 1 >, < V 2 >) has then a conservative form 
and expresses the balance of the chosen conservation laws. 

Remark 3.1. 

Let us come back now to Fig. [2] representing the spectrum of the operator —Q[x,T\/k. 
Let us consider the limit k — > now and consider the spectrum of —Q[x,T\/k on the space 
of 27r-periodic functions in this limit. We can say first af all that the sizes of the energy 
bands [E ,Ei], [E 2 ,E 3 ], [E^E^] tend to zero in this situation giving the "splitting" of 
the three localized quantum states of the corresponding decreasing (one-soliton) potential 
arising at k — > ojf) The potential <f>(kx; k, A, n) represents in this case a "lattice" of distant 
one-soliton solutions with the period ~ k^ 1 with k — > 0. 

Let us note now that the eigen- values of the operator —Q{ X ,T]/k for the energies E > E 6 
are double-degenerated on the space of periodic functions and represent the "boundaries 
of the gaps of zero width" in the spectrum of —Q[x,T\/k. It's not difficult to see also that 

3 The limit k — > of a one-phase solution of KdV gives a one-soliton solution corresponding to a 
reflectionless potential with one localized quantum state for the Lax operator (12.171) . The same solution 
gives a potential with three bounded states (one with E = 0) for the operator Q[x,t] given by p. 21) in 
the same limit. 
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the distance between these eigen-values decreases (~ k 2 ) in the limit k — > 0. Moreover, 
the size of the gap [E 5 ,Eq] (as well as the energy band [E^E^]) decreases also in this 
situation. As a result, we can see that another instability arises in our scheme for the case 
of the small k due to the large number of the "small" eigen-values of Q[x,t] m this limit. 
This fact means most probably that the averaging methods are not very applicable in the 
limit k — > where the "multi-soliton" description seems to give more adequate picture. 

At last let us say that we believe that it's enough to keep just the first (e 2 ) dispersive 
terms in system (I3.16P for the description of many oscillating regimes arising in the KdV 
theory. Finally we arrive at the system 

St = u{k,A,n) + e 2 u (2 ){k, A,n, ...) (3.22) 
k T = (u(k,A,ri) + e 2 u (2 )(k,A,n, ■ ■ ■)) x 

A T = a( 1) (k,A,n,kx,A x ,n x ) + e 2 a( 3 )( y k,A,n,...) (3.23) 
n T = ri il) (k,A,n,kx,Ax,nx) + £ 2 ri(3)(k,A,n,...) 

since we believe that it demonstrates already many essential features of the full system 

System (I3.23P should be considered as a system of differential equations in the ordinary 
sense, in particular, all the solutions of f !3.23|) are supposed to be well defined functions 
of X and T with some concrete behavior depending on the regime under investigation. 

In the next chapters we are going to consider the questions connected with the Hamil- 
tonian structures of system ( 13. 16f) which is the main subject of this paper. So, we will 
consider now the initial KdV equation as a part of an integrable hierarchy having two local 
Hamiltonian structures and discuss a possibility of the "averaging" of the Hamiltonian 
structures to obtain the Hamiltonian structures of the Dubrovin - Zhang type for system 

(GDI). 

4 The commuting flows and the Hamiltonian struc- 
tures. 

It is well known that the KdV equation represents the first nontrivial flow of the integrable 
KdV hierarchy generated by the higher KdV integrals 

/+00 
V v ((p,(p x ,...)dx 
-oo 

with respect to the Gardner - Zakharov - Faddeev bracket 

{<p(x),(p(y)} = 5'{x-y) 

or the Magri bracket 
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{<p(x),ip(y)} = 8"'{x-y) + - <p(x) 5'{x - y) + -ip x 8{x-y) 
All the higher KdV flows have the similar form 

ft- = f"(<P,<Px,<Pxx, ■ ■ ■ ) (4.1) 

and give an infinite set of commuting integrable flows. 

The commuting flows (14. If) can be also written in the "small dispersion" form 

ei p T „ = f v ((p ie <fix,e 2 ip X x,---) (4.2) 

which gives the commuting flows for the KdV equation written in the form (I3.18p . 

It is natural to expect then that the higher flows ( 14. 2 p of the KdV hierarchy gener- 
ate the commuting flows for the deformed Whitham system ( 13 . 16[) such that we get an 
"integrable" hierarchy starting from system (I3.16P on the "averaged" level. 

We have to introduce now the "extended functional space" M. = {ip(9, X)} consisting 
of smooth functions ip(8,X) which are 27r-periodic in 9 at every X. For our further 
purposes we need to introduce also a " submanifold" fC £ M. corresponding to the set 
of solutions ( I3.17P which will play the basic role in our considerations. Let us note here 
that all our considerations will be connected with the formal asymptotic series in the 
derivatives of parameters of one-phase solutions of KdV so we define also the submanifold 
fC in the same form, 1.6. clS 8b formal submanifold having the asymptotic sense. 

Thus, we define the submanifold K, in the space of functions <p(9,X) by the following 
rule: 

1) The function <p(9, X) belongs to the family /C if it represents one of solutions (I3.17p . 

i.e. 



<p(e,x) 



$ + 9, k, 4 n) + E e ' $ « {^T 1 + 9 > ^ A > n l ' X 



with some functions (S(X), A(X),n(X)) where k(X) = Sx(X); 
2) We put the following relation between the functions S(X) and 

S(X) = - / sgn(X -Y)k(Y)dY (4.3) 

The functions (k(X), A(X),n(X)) play the role of "coordinates" on the submanifold 
/C, so we consider K, a manifolds parametrized by three functional parameters. 

Let us formulate here the Theorem which connects the higher flows (14.21) with the 
commuting flows of system ( 13.161) . 



4 Let us assume here that the relations k{X) — > for X — > ±oo are imposed. However, the procedure 
will give us a local deformed Poisson bracket on the space (k(X), A(X),n(X)), so this condition will not 
be important in fact for the final result. 
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Theorem 4.1. 

Every higher KdV flow (f^.ffi ) leaves invariant the family of formal solutions j3.11\) 



and generates a commuting flow for the deformed Whitham system l}3.16}) which can be 
represented in the same graded form 

S TV = u} v (k,A,n) +J2e 2l u u m (k,A 1 n 1 k x ,A x ,n x ,...) (4.4) 

i>i 



k T v = [u)"(k,A,n) + y^^ 21 U( 2l] (k, A, n, k x ,A x 



i>i 



A T „ = ^e 2/ afa +1) (k,A,n,k x ,A x ,n x ,. . .) (4.5) 



l>0 

n TU = ^2e 2l r] u (2l+1) (k,A } n,k X} A X} n x , 

l>0 



as system 113.16}) . 
Proof. 

Let us consider the formal asymptotic series 



(4.6) 



where every function ^(i)(9, X, T, T v ) is a local functional of k Q (X,T) = So X (X,T), 
Aq(X,T), no(X,T) and their X-derivatives which is polynomial in derivatives and has 
degree I according to the gradation rule we introduced above. We require that series ( 14. 6 j) 
coincides with asymptotic series ( 13. 17ft with the same parameters ko(X,T), A (X,T), 
n (X,T) for T v = 

tf (O )(0,X,T,O) = H9,k ,A , n ) , tf (o (0, X,T,0) = $ {l) (6,X,T) 

and satisfies the higher KdV equation (14.21) for T u > 0. 

After the substitution of (14.61) into (14.21) in the graded form we get a chain of evolution 
equations on the functions ^(i)(0, X, T, T u ) at every degree I 

V (l) (e,X,T,T u ) = A {l) (y,y e: ty x ,...,ko,Ao,n ,kox,A ox ,nox,---) (4.7) 



dT 

where every A(/) depends only on ^i s ) with s < I. 

It's not difficult to check the following relations for T v = 



* (O (-0,X,T,O) = (-l) l V [t )(e,X,T,0) , A (l) (-9,X,T,0) = (-l) m A (O (0, X, T, 0) 

(4.t 
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I > 0. 

Let us assume for simplicity that all the systems (14. 7p have smooth solutions on some 
interval T v G [0, 5} with our initial data such that we get a unique formal series (14. 6[) 
satisfying our requirements on the same interval. Since equation (14. 2 p gives a commuting 
flow for the KdV equation (13.181) we get that series (14. 6p gives a formal solution of ( 13. 18j) 
at every T u G [0,5]. However, series (I4.6p can not be considered as the asymptotic series 
(I3.17P for T v > since the normalization conditions (12 .25 p . (I3.12p will be in general 
destroyed by the evolution systems (l4.7p FI 

Nonetheless, series (14 .6p can be represented in form ( 13.17P after a redefinition of pa- 
rameters 



(k (X,T),A (X,T),n (X,T)) -> (k(X, T, T"), A(X, T, T"), n(X, T, T")) (4.9) 

and a re-expansion of (14.61) in the new graded form. It is convenient then to represent 
the redefinition of (k, A, n) in the differential graded form (14.51) which gives the required 
evolution system on family (I3.17p . 

Let us discuss finally the possibility of a construction of the required system (14.51) on 
the space of parameters (k(X), A(X),n(X)) which will prove the Theorem. Indeed, the 
function n(X,T,T u ) is given by the integral 

r2n in 

/ <p(e,X,T,T») - 

Jo 271 
according to the definition, so we get immediately the graded equation 

1=0 J{i 

which makes satisfied the second relation (I2.25P for T u > 0. 

Let us consider now the first relation (12.251) and relation (13.121) . We have to find now 
two more functions S(X, T, T u ), A(X, T, T u ) such that the function $(6*, Sx, A, n) satisfies 
the conditions 



r-2n 



'0 



^ e ' i>o ' 71 



fS(X J T 2 T± + x 



x 



5 The function ^( ) {^i Xi T, T v ) will remain the one-phase solution for T" > in this situation, however, 
the normalization ^( )e(0, X, T, T v ) — will be also destroyed by the higher KdV flow on the one-phase 
solutions. 



25 



to be an appropriate main term in the "re-expanded" series H4.6H . 

We assume also S(X,T,0) = S (X,T), k(X,T,0) = k (X,T), A(X,T,0) = A (X,T), 
n(X, T, 0) = uq(X, T) according to our scheme. 

Differentiating the first relation with respect to T v at T v = we get 

( - S T » $ee + k T , $ 0fc + Ap, <5> eA + n T » $ e „ J ^ e l * w — + 

/27T JO 
$^e J A (0 | (T , =0 )— = (4.10) 

We are going to obtain a graded linear system for the determination of the time 
derivatives St», A^ in the graded form. Using the facts 

tt (o) (0,X,T,O) = $( g(I ' T) +M,A, B 



e 

A w (»,Jf,r,0) = * 9 + A, » 

where a; y (fc,v4,n) is the frequency corresponding to the flow f u on the space of the one- 
phase solutions of KdV we get from equation (14.101) St« = u) v (k,A,n) at T u = in the 
main approximation. 

Using also relations (14. 8p we can write actually 

S T u = u"{k,A,n) + £>(e 2 ) 

at T u = for the derivative St»- 

After the differentiation of the second relation w.r.t. T u at T v = we get the following 
relation 



/*2tt ^ ;-27r / j \ _^ 

At* / ii<& A — = I (- S T » £io + k T » £ifc + ^t- 6a + ny* £m ) XI e ' $ W o - 





A(2{+1)I(T»=0) - ^T"$fc - ™T"$n ^~ (4.11) 



at T 1 ' = 0. 

The function J Q n ^ $^4 d8/2n is a strictly positive function on the space of parameters 
(k,A,n). Using this fact it is not difficult to see then that the form of the linear system 
(I4.10p - (l4.1ip defines the unique representation of the derivatives St», A^ at T v = in the 
graded form being uniquely resolvable at every step of the determination of Sjl and Aji . 
Using also relations (14. 8 p it's not difficult to prove that we obtain the "purely dispersive" 
system (14 .5p in this situation. 
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For < T u < 5 system (I4.10p - fl4.lip is still resolvable with respect to the derivatives 
St", At», such that we have 





= t% W (20( T "' n °> k ° x > A ox,n ox , ■■■ ] 












= ( ^2 g2/ w (2o(^> fc °> n °' A ° x ' n °^' • • 


i 




\Z>0 




Aj>v 


= ^ a (22+l)(^' ^0) A); ^0; ^0X, ^0X, • 


••) 




l>0 





n T „ = J ^2e 2l r]'( 2l+1) (T u ,k ,A ,n , k ox ,A 0X ,n ox ,...). 

where all the functions oj" 2 i)i a (2Z+i)' V(2i+i) become dependent on T v and do not coincide 
with the functions from ( 14. 5 p since all the functions ^n\{9, X, T, T v ) become different 
from <&(j')(0,X,T). However, if we represent the solutions of this system, say, in the 
formal graded form we will be able to "re-expand" formal solution (14. 6p according to 
change of parameters ( 14. 9p . 

Thus we can write now the new formal graded expansion for series (14. 6 p 

<p (9, X, T, T") = J2 e ' *(0 + 9, X, T, tA (4.12) 

i>o ^ e ' 

according to change of parameters of expansion (14.91) at every T v . System (14. 7p can also 
be easily rewritten for the functions ^(i)(9, X, T, T u ) 



—V^X^T) = A {l) (T\V,*e,*x,...,k,A,n,k x ,A x ,n x ,...} 

using system f)4.12p in this situation. As a result we will get asymptotic series f )4.12p 
satisfying all the conditions (I')- (III'). 

Finally we get that the asymptotic series (14.121) gives a formal solution of (13 . 18[) satis- 
fying all the conditions (I')-(III') at < T u < 5. As we saw above, solutions (14. 12j) should 
coincide in this case with the formal graded solution (I3.17P so we get the invariance of the 
family (I3.17P under the higher KdV flows. The evolution of parameters (k, A, n) is ruled 
then by system (14. 5p for all T v > and the commutativity of (14. 5p with (I3.16P follows 
directly from the commutativity of (I4.2p and (I3.18p . 

At last, let us note now that for our conditions k — > 0, X — > ±oo we have also 
St"{X) — > 0, X — >■ ±oo, which gives also the conservation of condition (14. 3p in our 
situation. 

Theorem is proved. 
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It's not difficult to see also that system ( 14. 5 j) coincides with the deformed Whitham 
system for the higher KdV flow (14. 2 p defined by the same normalization conditions (I')- 
(III'). 

Let us discuss now the Hamiltonian properties of system (I3.16P following from the 
Hamiltonian properties of the KdV equation (13.1 8ft . According to the general ideology 
of the deformation of systems of Hydrodynamic Type we will assume the existence of 
Hamiltonian structures for the deformed Whitham system given by the deformations of 
the Hamiltonian structures of Hydrodynamic Type, i.e. the Hamiltonian structures having 
the form 

{U V {X),U»{Y)} = {ir(X),U"(Y)} + 

k 

+ E fi (t)s( U > U *>'--> U (^) (4.13) 

k>2 s=0 

where all B v ,ji s are polynomial w.r.t. derivatives XJx, • • • , V(k-s)x and have degree (k — s). 

We call deformations of Hamiltonian structure of form (14.131) the deformations of 
Dubrovin-Zhang type. Bracket (14. 13|) gives a deformation of the local homogeneous 
bracket of Hydrodynamic Type (Dubrovin - Novikov bracket) which according to the 
definition has the following form 

{U U {X),U' M {Y)} = g v "(U)8'(X -Y) + b u ^(V)U^5(X -Y) (4.14) 
The corresponding Hamiltonian operator J UfJ - can be written as 

ju» = ^ (U) jL +6 ^ (U)f/ A 

and is homogeneous w.r.t. transformation X — > aX. 

Every functional H of Hydrodynamic Type, i.e. the functional having the form 



H = / h(U) dX 

J — oo 

generates the system of Hydrodynamic Type 

XJ V T = V»(U)U£ , = 1,...,JV (4.15) 

where V"(U) is some N x N matrix depending on the variables U l , . . . , U N according to 
the formula 

The DN-bracket ffl~T4j) is called non-degenerate if det ||^(U)|| ^ 0. 
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As was shown by B.A. Dubrovin and S.P. Novikov the theory of DN-brackets is closely 
connected with Riemannian geometry ([TTJ [121 m fact, it follows from the skew- 

symmetry of (I4.14p that the coefficients g Ufl (XJ) give in the non-degenerate case the con- 
travariant pseudo- Riemannian metric on the manifold Ai N with coordinates (U 1 , . . . , U ) 
while the functions r^ A (U) = — p MQ ,(U) 6°"(U) (where g utl (U) is the corresponding met- 
ric with lower indices) give the connection coefficients compatible with metric g Ufl (XJ). 
The validity of Jacobi identity requires then that g Ufl (XJ) is actually a flat metric on the 
manifold Ai N and the functions r^ A (U) give a symmetric (Levi-Civita) connection on 

m n ([mm us]). 

In the flat coordinates n 1 (U), . . . , n N (U) the non-degenerate DN-bracket can be writ- 
ten in constant form: 



K(X)X(y)} = e u 5^5'(X-Y) 



where e u = ±1. 
The functionals 



/+oo 
n u (X)dX 
-oo 

are the annihilators of the bracket (14.141) and the functional 



1 r+°° N 
P = ~ £V (n*(X)) 2 dX 

Z J -°° u=l 



is the momentum functional generating the system U!f = U\ according to f)4.16p . 

The Symplectic Structure corresponding to non-degenerate DN-bracket has the weakly 
nonlocal form and can be written as 



Q vlt (X,Y) = e v 5 vll v{X-Y) 
in coordinates n u or, more generally, 

^(x,^ = Xy££(*MJr-y)|£(y) 

A=l 

in arbitrary coordinates U u . 

Let us mention also that the degenerate brackets f)4.14p are more complicated but also 
have a nice differential geometric structure ([33]). 

Brackets f)4.14p are closely connected with the integration theory of systems of Hy- 
drodynamic Type f)4.15p . Namely, according to conjecture of S.P. Novikov, all the di- 
agonalizable systems (I4.15P which are Hamiltonian with respect to DN-brackets (14.141) 
(with Hamiltonian function of Hydrodynamic Type) are completely integrable. This con- 
jecture was proved by S.P. Tsarev ([55]) who proposed a general procedure ("generalized 
Hodograph method") of integration of Hamiltonian diagonalizable systems (I4.15p . 
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In fact Tsarev's "generalized Hodograph method" permits to integrate the wider class 
of diagonalizable systems (I4.15P (semi-Hamiltonian systems, |65j) which appeared to be 
Hamiltonian in more general (weakly nonlocal) Hamiltonian formalism. 

The corresponding Poisson brackets (Mokhov - Ferapontov bracket and Ferapontov 
bracket) are the weakly nonlocal generalizations of DN-bracket (I4.14p and are connected 
with geometry of submanifolds in pseudo-Euclidean spaces. Let us describe here the 
corresponding structures. 

The Mokhov - Ferapontov bracket (MF-bracket) has the form ([57]) 

{U v (X),U fl (Y)} = g^(XJ)5\X-Y) + b^\\J)U x 5(X-Y) + cU x v(X-Y)U$ (4.17) 

As was proved in [57] the expression (14.1 7p with det \ \g uti (\J)\\ ^ gives the Poisson 
bracket on the space U U (X) if and only if: 

1) The tensor g u,x (JJ) represents the pseudo-Riemannian contravariant metric of con- 
stant curvature c on the manifold Ai N , i.e. 

2) The functions T^ A (U) = —g^TJ) b^"(U) represent the Levi-Civita connection of 
metric ^(U). 

The Ferapontov bracket (F-bracket) is more general weakly nonlocal generalization of 
DN-bracket having the form ( [251 1261 [271 [28] ) : 

{U u {X),U fi (Y)} = g ulx (U)6'{X -Y) + 6^(U) U x 6{X - Y) + 

+ J>^ fe)A (U) U x v{X - Y) ^ )5 (U) U S Y (4.18) 
k=i 

e k = ±1, v,/jl = 1, . . . , N . 

The expression (14.181) (with det \ \g v ^(U)\\ ^ 0) gives the Poisson bracket on the space 
U V {X) if and only if [28]): 

1) Tensor g utl (XJ) represents the metric of the submanifold M. N C K N+9 with flat 
normal connection in the pseudo-Euclidean space K N+9 of dimension N + g; 

2) The functions r^ A (U) = —g fJia (XJ)b'^ u (XJ) represent the Levi-Civita connection of 
metric ^(U); 

3) The set of affinors {w;^ A (U)} represents the full set of Weingarten operators cor- 
responding to g linearly independent parallel vector fields in the normal bundle, such 
that: 

MUK fc)/l (u) = ffMr(UK*>(u) . v, w f fc)A (u) = v A ^(u) 
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^(U) = (^)a(U)^ )?? (U) - ^ )A (U)^(U)) 

k=l 

Besides that the set of affinors W(k) is commutative [itf(fc), = 0. 

As was shown in [26] the expression (I4.18P can be considered as the Dirac reduction of 
the Dubrovin-Novikov bracket connected with metric in K N+9 to the manifold Ai N with 
flat normal connection. Let us note also that MF-bracket can be considered of 
the F-bracket when Ai N is a (pseudo)-sphere S N C K N+1 in a pseudo-Euclidean space. 

The Symplectic Structures Q UfJj (X, Y) for both (non-degenerate) MF-bracket and F- 
bracket have also the weakly nonlocal form ([53]) and can be written in general coordinates 
U v as 

n vtl (X,Y) = £^_(x M x_y)_(y) 

s=l 

where e s = ±1 and the metric Gu in the space K N+9 has the form Gu = diag(ex, . . . , e^ +g ) 
The functions n l (\J), . . . ,n N+9 {\J) are the "Canonical forms" on the manifold Ai N and 
play the role of densities and annihilators of bracket (14. 18j) and " Canonical Hamiltonian 
functions" (see [53]) depending on the definition of phase space. In fact, the functions 
n s (XJ) are the restrictions of flat coordinates of metric Gu giving the DN-bracket in K N+9 
on manifold M N . The mapping M N -> E N+9 : 

(U\...,U N ) (n\U),...,n N + 9 (V)) 
gives locally the embedding of Ai N in K N+9 as a submanifold with flat normal connection. 

All the brackets ( 14.141) . ( 14.171) . ( 14.181) are connected with Tsarev method of integration 
of systems (14.151) . Namely, any diagonalizable system (I4.15P Hamiltonian w.r.t. the (non- 
degenerate) bracket (14.141) . (14.171) or (14.181) can be integrated by "generalized Hodograph 
method" . 

We will not describe here Tsarev method in details. However, let us point out that 
"generalized Hodograph method" and the HT Hamiltonian Structures were very useful 
for Whitham's systems obtained by the averaging of integrable PDE's ( [701 1291 [TT] |4~21 l4"3~t 

mm)- 

The Hamiltonian approach to the Whitham method was started by B.A. Dubrovin and 
S.P. Novikov in [11] (see also [121 US]) where the procedure of "averaging" of local field- 
theoretical Poisson bracket was proposed. The Dubrovin - Novikov procedure gives the 
DN-bracket for the Whitham system (14.151) in case when the initial system is Hamiltonian 
w.r.t. a local Poisson bracket 

V(v)} = Y. B U^^---)s {k \*-y) 

k>o 
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with local Hamiltonian functional 

/+oo 
h(<p,(p x ,...)dx 
-oo 

The method of Dubrovin and Novikov is based on the presence of iV (equal to the 
number of parameters IF of the family of m-phase solutions) local integrals 



F = J V v (cp,cp x ,...)dx 
commuting with the Hamiltonian function and with each other 



(4.19) 



{F,H} = , {F,F} = (4.20) 

and can be formulated in the following form: 

We calculate the pairwise Poisson brackets of the densities V v in the form 

{v»(x), v»( y )} = A 7iv, v>«, ■ ■ ■ )$ {k \x - y) 

where 

A^(<p,<p x ,...) = d x Q^(ip, ( p x ,...) 

according to ( 14.201) . Then the Dubrovin- Novikov bracket on the space of functions U(X) 
can be written in the form 

{IF(X), IF(Y)} = {AT){U) 5'{X -Y) + -^-1% S(X - Y) (4.21) 
where (...) means the averaging on the family of m-phase solutions given by the formula: 



and we choose the parameters U u such that they coincide with the values of F on the 
corresponding solutions 

ir = (p»(x)) 

This procedure was generalized in [52] for the weakly nonlocal Hamiltonian structures. 
In this case the procedure of construction of general F-bracket (or MF-bracket) for the 
Whitham system from the weakly non-local Poison bracket for initial system was pro- 
posed^] In [54] the procedure of averaging of the weakly-nonlocal Symplectic structures 
was also suggested. 

6 The final proof of the Jacobi identity for the bracket given by the Dubrovin - Novikov procedure was 
given in |51) . 
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Here we are going to consider the construction of bracket (14. 13j) for the deformed 
Whitham system (I3.16p . Being considered for the KdV case the procedure will have in 
fact general character and can be considered as a generalization of the Dubrovin - Novikov 
procedure for the deformed Whitham system in general case. 

We are going to use the Dirac restriction of a Poisson bracket on a submanifold to 
establish the procedure of the construction of a Poisson bracket of form ( I4.13P for the 
deformed Whitham system (13 . 1 6[) which we call the "averaging" of a Poisson bracket for 
the deformed Whitham systems. Let us first introduce the Poisson brackets 

{<p(6,X),<p(9',Y)} = e6(9-9')5'(x-y) (4.22) 

and 

{<p(6,X),<p(e',Y)} = (4.23) 

= e 3 5(6 - 6') 5"'(X - Y) + e - ip(X) 5(6 - 6') 5'(X — Y) + e - tp x 5(6 - 6') 5(X - Y) 

3 3 

which correspond to the Gardner - Zakharov - Faddeev and the Magri brackets on the 

extended phase space (p(6, X) periodic in 6 with the period 2ir. Easy to see that both the 

expressions ( I4.22[) and (I4.23P give Poisson brackets on the extended functional space. 

We have to consider now the " subspace" in the extended functional space correspond- 
ing to the full family ( I3.17P parametrized by three functional parameters (S(X), A(X),n(X)). 
We will call now the Dirac restriction of bracket (I4.22p or (I4.23P on the submanifold 
corresponding to the full family of solutions (I3.17P the averaging of the Gardner - Za- 
kharov - Faddeev bracket or the Magri bracket giving a Poisson bracket for the deformed 
Whitham system (I3.16p . 

Let us remind that the Dirac restriction of a Poisson bracket on a submanifold Af k C 
M. n is connected with the special choice of coordinates in the vicinity of the subman- 
ifold M k which are divided to the "coordinates on the submanifold" (U 1 ,...,U k ) and 
the constraints (g , . . . ,g n ) which define the submanifold J\f . It is assumed that the 
submanifold M k is given by the conditions 

g l (x) = , i = l,...,n — k 

while the k functions f/ 1 (x), . . . , U k (~x) on M. n play the role of coordinate system on M k 
after the restriction on this submanifold. 

If the Hamiltonian flows generated by the functions f/- 7 (x) leave the submanifold M k 
invariant, i.e. we have 

{fP(x),(f(x)} = for g(x) = 

then the pairwise Poisson brackets of functions t/ J (x) give a Poisson tensor after the 
restriction on M k with coordinates (U 1 , . . . , U k ) which is called the Dirac restriction of 
the Poisson bracket {...,...} defined on Ai n on the submanifold M k C Ai n . 

7 Let us note that all the constructions are considered here just on the level of the formal asymptotic 
series. 
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In general, according to Dirac procedure, if we have some constraints g*(x) which 
define a submanifold M k and some functions (x) giving a coordinate system on Af k we 
have to find k linear combinations f3 J s (XJ) g s (x.) at every point of M k such that we have 
for the functions 

IP(x) = *P(x) +#(U)«7'(x) , j = l,...,* 

the relations 

{^(x),^(x)} = at g(x) = 

The functions f/ J (x) have the same values as the functions £7 J (x) at the points of M k 
and we can then define the Dirac bracket {...,... }d on M k by the formula 

{U\W} D = {^(x),^(x)}U(U) 
The functions /3s (U) are defined from the linear system 

tf(x),g'(x)}\j< r .Pi{V) + {/(x),^(x)}|^ = , i = l,...,n-k 
and we can also write 

{U\W} D = {^(x),^)}^ - /3:(U){/(x),^(x)}U fc ^(U) 

for the Dirac bracket on A/" fc . 

Let us describe now the Dirac procedure in our situation. 

First, let us introduce new coordinates on the submanifold fC corresponding to solu- 
tions (13.171) based on the conservation laws of the KdV equation (I3.18p . 

Let us choose three integrals of the KdV equation such that their values on the family of 
one-phase solutions of KdV are functionally independent. In our case it is most convenient 
to take the integrals 

f +QO v 2 f +co fv 2 <{%\ 

yds, h = -d X) I 2 = j ^ [ Y - Y )dx 

The integrals transform naturally to the integrals of the KdV equation (I3.18P on the 
extended phase space 

+ °° dX— I f + °° !^ ^ dX — 

oo Jo ^ 2vr ' 1 J 2 2vr 



J -oo JO 

Let us introduce now the functionals 

r2w 



2 rA 2vr 



/ df) 
V v {ip,eip X) ...) — , i/ = 0,l,... 
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i.e. 

u"(x) = f\(e,x) g , u\x) = \j\\e,x) d ± 

^w = jf (^ ,jr) -^^)i (424) 

and consider the values of the functionals U U (X) on the submanifold /C. 

It is easy to see that the values of U U (X) on /C are equal in the main approximation 
to the values of the functionals l v on the one-phase solutions of KdV with the parameters 
(k(X), A(X) ) n(X)) and have in general higher corrections polynomial in derivatives of the 
functions k(X), A(X), and n(X). It's not difficult to see also that the higher corrections 
to U"(X) contain only even degrees in the expansion w.r.t. the derivatives kx, Ax, fix, 
. . . for our choice of the initial phase of the functions k, A, n) view the statements of 
the Lemma 3.1. Thus for the functionals U°(X), U 1 (X), U 2 (X) we can write 

U\X) = (MX) ee n(X) 

U\X) = (v, 2 ) (X) + J^Ul (4-25) 

S>1 

S>1 

where 

/27T 
$>(6,k,A,n) — = n 

/■27T in 

<^}„ = J <t 2 (e,k,A,n) 

and the values U\ s , U$ s are graded polynomials in the derivatives of k, A, and n having 
degree 2s. 

Since the values of 1°, J 1 , I 2 are functionally independent on the space of one-phase 
solutions of KdV we can write the "inverted series" for the functions k(X), A(X), n(X). 
We have to change the gradation rules now such that we will define the gradation degree 
with respect to the X-derivatives of the parameters (U°, U 1 , U 2 ) instead of (k,A,n). So 
we can write now 

k(X) = ko(U°(X),U\X),U 2 (X)) +J2k2s)([U°,U\U%X) 

S>1 

A(X) = A (U°(X),U\X),U 2 (X)) + J2 A (2s)([U°,U\U 2 ],X) (4.26) 

S>1 
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n{X) = U°(X) 

where ko(U°, U 1 , U 2 ), A (U°, U 1 , U 2 ) are the exact "one-phase" expressions for the param- 
eters k and A in terms of (U°, U l , U 2 ) and the functions k( 2s ), ^.(2s) are graded polynomials 
in the derivatives of (U°, U 1 , U 2 ) having degree 2s. 

Using relations (I4.26j) we can re-expand also solutions (13 . 1 Tl) as graded series with 
respect to the X-derivatives of the values of the functionals U°(X), U 1 (X), U 2 (X) on K. 
at every time, such that we have 



<j>(0,X,T) = $° (^^l + 9,U ,U\U^j + J^e l ^ l) (^^ + e,X,T S j (4.27) 

where all $^ are graded polynomials of (U x , U x , U x , . . . ) of degree I. 

The function $ U (6 I , U°, U 1 , U 2 ) represents the exact one-phase solution of KdV de- 
pending on parameters (U°, U 1 , U 2 ) and we have by definition 

<t> v {e,u ,u\u 2 ) = $(e,k(u),A(u),n(u)) 

According to our approach we will assume that series (14.271) and (13.171) are equiv- 
alent representations of formal asymptotic solutions ( 13.1 7ft connected by change of the 
asymptotic functional parameters ( I4.25p - (l4.26p . Let us note also that due to the form of 
relations (I4.25[) - fl4.26[) the symmetric properties of the functions $^(6*, X, T) remain the 
same as for the terms of series ( I3.17p . i.e. we have 



^ 2s) (-9,X,T) = ^ 2s) (9,X,T) , <S>? 2s+1) (-9,X,T) = -$f 2s+1) (9,X,T) , s>0 

We can assume in the same way that the functions of the submanifold K. are represented 
now by the asymptotic series 



tp{0,X) = $ u (^1 + 9,U°,U\U 2 ) + Y j e l ^ ) (^ + 9,[U°,U\U 2 ],x) 
^ e ' i>i ^ e ' 

so the functionals U°(X), U 2 (X) play the role of coordinates on this submanifold 

and these are exactly the functionals we are going to use for the Dirac procedure. 

Let us introduce now the system of " constraints" which defines our submanifold K. in 
the functional space. For our purposes it will be convenient to write the constraints in 
the following form: 

Let us denote 



i>(9,[U],X) = $ V (9,U,X) + ^e<$g(0,[U],X) 

1>1 
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where the notations U(X) = (U°(X),U 1 (X),U 2 (X)) denote now the functional fTC24l 
defined on the full functional space {(p(9,X)}. We introduce now the constraints g(9,X) 
by the formula 

g(9,X) = <p(9,X) -J g[U | (I) +g,[U],lj (4.28) 
as the functionals on the space M.. Easy to see that the relations 

9(0, X) = 

define then exactly the " sub- manifold" /C we consider here. 

However, the set of constraints ( I4.28P is certainly not independent in the ordinary 
sense. Namely, in the full analogy with the finite-dimensional case the following relations 
take place identically for the "gradients" 8g{9 , X) /Sip(9' ,Y) on the "sub- manifold" /C: 

Nevertheless, it will be convenient for us not to choose an independent system of 
constraints and to keep constraints (I4.28P for our purposes, so we have to remember the 
presence of relations (I4.29[) for system (14.281) . 

For the Dirac restriction of bracket (I4.22p or (I4.23P on the submanifold /C we have 
to modify now the functionals U°(X), U 1 (X), U 2 (X) by the linear combinations of con- 
straints g(9, X) 

U»(X) = IT(X) + | + ° °jT g(9,Y) 0-f^Ml + 9,[U],Y,x'\ ^-dY 

such that the functionals U V (X) leave invariant the submanifold K. in the corresponding 
Hamiltonian structure and then to use the functionals U U (X) for the construction of the 
Dirac bracket on /C. The functions (3"(S(Y)/e + 9,Y, X) should satisfy the relation 

J + °° j\g(9,X),g(9>,Z)}f3» f^p- + 9',Z,Y^j ^-dZ + {g(9,X),U»(Y)} = 

(4.30) 

on /C and are defined at every "point" of K. modulo the linear combinations of the functions 
5U' 1 (W)/5tp(9, Y) view the original dependence of constraints ( I4.28p . 
The Dirac bracket on the manifold K, can be defined by the formula 

{U»(X),U»(Y)} D = {U U (X), U ti (Y)}\) C — (4.31) 
- J F + 9, Z, {g(9, Z),g{9', W)}\ K fi» + 9\ W, Y^j ~dZdW 

so the procedure gives a unique definition of the bracket {U U (X), U^{Y)} D . 
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To get a local deformed Poisson bracket on K. we will try to find the functions 
f3 u (9,Y,X) in the form 



P v (0,Y,X) = Y,z S PU e > Y i X ) ( 4 - 32 ) 



S>1 



where the functions l3Y a \(8, Y, X) are represented as the local distributions 

s 

v {s) (9,Y,X) = ^p^ifi^S^Y-X) (4.33) 

p=0 

having gradation s assuming that the derivatives of the delta-function 5^ (Y — X) have 
degree p by definition. 

Thus, we assume that all the functions 0^ p (9,Y) on K, are local functionals of 
(U°(X),U 1 (X),U 2 (X),U X ,U X ,U X ,...) at every 0, polynomial in derivatives 
(U X ,U X ,U X , . . .) and having degree s — p according to our previous definition. This 
structure of fl u (9,Y, X) is obviously equivalent to the statement that the functionals 

x» 

U U {X) q(X) dX 

with a "slow" function of X q(X) can be modified with the aid of a linear combination 
of constraints (14.28)1 with the coefficients 



dXP 

S >1 p=0 

to leave the submanifold K, invariant. According to this scheme the derivatives d s q/dX s of 
the slow function q(X) have degree s as well as the derivatives of the parameters U°(X), 
U\X), U 2 (X). 

Finally, we have to study now system ( I4.30P for the cases of the Gardner - Zakharov - 
Faddeev bracket and the Magri bracket to investigate the possibility to find the functions 
(3 U (8,Y, X) in form (I4.32l) - (l4.33p . Let us formulate here the following Theorem. 

Theorem 4.2. 

Both for the Gardner - Zakharov - Faddeev bracket and the Magri bracket for KdV the 
functions P u (0, Y, X) can be found in form \J^.32 )- ^.3S ) on the family /C. Thus, the Dirac 



restriction of both the brackets on the family K, has the local deformed Hydrodynamic form 
Iji4-13\ ) which gives two deformed Hydrodynamic Type brackets for the deformed Whitham 
system A3. 16}) . 

Proof. 

Let us analyze equations (I4.30P for the case of the Gardner - Zakharov - Faddeev 
bracket and the Magri bracket. We have first on the family /C 
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{g(e,X),g(9',Z)}\ lc = {<p(9, X), p(6', Z)}\ 



oo </ — oo 



summation over repeated indices) Jf] 
In the same way 

{g{0,x),v(Y))\ K = {p(0,x),cr(Y)}| K - 

Let us say now some words about the Poisson bracket {k(W),U u (Y)}\ic- As we saw 
already the functional 



/+oo 
U U {Y) dY 
■oo 



leaves invariant the submanifold /C so the Poisson bracket {k(W), I u }\ic should give exactly 
the Whitham evolution of the functional fc([U], W) corresponding to the z/-flow of the KdV 
hierarchy. So we have 

{k(W),n\K = e[uj"(k,A,n) + Y j e 2l uj'( 2l) (k,A,n,...) 



eL"(k,A,n) + J2 

V i>i 



w 



with some functionals uY^JJJ, Ux, • • • ) according to (14. 5p . 

According to the structure of the bracket {k(W), U v {Y)}\k. we should have then 



{k(W),U»(Y)}\ K = 6 w£(U(W)) + ^e 2 '^ 20 (U,U W) ...) 8(W-Y) + 

V i>i J w 

S>1 

8 Let us note that we assume the differentiation S/SU in the sense of the values of functionals U U (X) 
on the family K. while we treat U V {X) inside the brackets as a functional on the whole functional space. 
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where k"^ are some local functionals of (U, Uw, • • • ) given by sums of terms of degree 
> 0. 



We can write then 



1 

2e 



+00 

dWMO, [U],X) sgn(X- WO {k(W),U"(Y)}\ K 



= M0,M,X) (u> v (V(X)) + |> 2 ^ 20 (U,U x ,...)j S(X-Y)- 
l - J + °°dW MO, [U], X) sgn{X-W) (uZ(U{W)) + £ e 21 ^(U, V w , . . . ) j S'(W-Y)+ 
I l +O °dWM0,M,X)sgn(X-W) (^e s K^(V,V w ,...)S^(W-Y)j 



+ 2e 



In the same way we obtain 



M^x),r}k = ^(^p- + 9,[v],x^j ^(u(x)) + |J e 2 '^ 20 (u,u x ,...)j + 

+ J + °° Vt/M + 0, [U] , X, W^j {U»(W) ,r}\ K 

where 

So, from the structure of the bracket {(/?(#, X), C/^(y)}|^ we can conclude 

xl(^p- + o,x,Y^j = Me,x),u^Y)}\ K = 
= i>e + e, [u] , x^j (ut(v(x)) + 1] e 2/ afo (u, Ux, . . . ) j 5(x - y) + 

+ J^dWiPu, (^l + e,[V},X,W^j {U»{W),U V (Y)}\ K + (4.34) 

for some local functionals A^(#, [U], X) on /C, polynomial in the derivatives (Ux, Uxx, • • • ) 
and given by sums of terms of degree > 0. 
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In the same way we put 



\ [YX.^-9. ) = {U»(Y).f((>.X)}U 



= -S(Y- X) (u%(V(X)) + £ e 21 u>fo(U, U x , . . . ) j ^ + 0, [U], + 

+ dW {U V {Y), U»{W)}\ K ifo, + 9, [U], X, + (4.35) 
+ £ e * (y - X) A# + 61, [U] , x) 

s>l ^ 6 ' 



s 

Let us denote also 



C L (^p- + e,X,Y,' S ) = M9,X),k(Y)}\,c 
C R ( Y ,X,^l + 9,) = {k(Y),ip(9,X)}\ K 



We have now 



or(f!W + e,x,Y) = {<j(0.x).i »(Y)}\k = 



(^ + o,x,y) 

= {ip{9,X),U»{Y)}\ K - j +0 °dW^ (^l + 9,[V],X,W^j {U»{W),U»{Y)}\ K - 
~t J + ^ dW ^o(^p- + e,[V],X^ sgn(X-W) {k(W),U»(Y)}\ K = 

= E eS A w (—+0, [u]. x ) 5(S) (^ - n + 

s>l ^ e ' 

+ l - J + °°dW^ (^l + e,[XJ],X^ sgn(X-W) (u» (W) + £ e 2l u\ 2l) {W) j 

dWfr(?^ + 0,\V\,x) sgn(X-W) ^2e s K^(W)5^(W-Y) 

Using the same arguments we obtain that for the case of the Gardner - Zakharov 
Faddeev bracket we have the following equation for the functions j3"(0,Z,Y): 



/+oo i>2tt rn/ 
J L(9,9',X,Z) (3»(9',Z,Y) — dZ = a»(9,X,Y) (4.36) 
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where 

L(e,e',X,Z) = k5'(6 -6')5(X - Z) + e5(e-e')8'(X -Z)- 

1 /" + oo 

"27 7 ^CL(#,X,^)sgn(Z-^)^',[UU)- 



1 

27 

-oo />+oo 



oo 
+oo 

dWMO, [U],X) sgn(X- ^) Cr(W,Z,9') + 



+ £ / J LI dWdV 6 ^uhw) ] mwl k{m,c sgn {z ~ v) M *> [UU) + 

+ ^ J + _™ dW dV MO, [U], X) sgn (X - ^ 7 (^)}k ^§j^ + 

^ J— oo J ~oo 

Let us note now that the bracket {[/"(X), has the order 0(e) and, besides 

that, it's main term in the e-expansion coincides precisely with the Dubrovin - Novikov 
bracket defined above. 

Let us put now the additional condition 

J M0,[U],Z) P V (6,Z,Y) — = (4.37) 

which will be confirmed aposteriori for our /3 U (9, Z, Y). We can reduce then the operator 
L(9,6',X, Z) to the form 

L eff (e,e',x,z) = k5\e -e')5(x - z) + e8(e -o')5\x - z) - 

i /-+oo 

__ y dwM0,\u\,x) sgn(x-w) uw,z,e') + 



+ - / / dWdy^(0,[U],X)sgn(X-W) {A;(W),^(y)}|x; 

^ e </-oo ./-oo 
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Let us define now the functions (3 U (9, Z, Y) as the solutions of the equations 

"+oo /-27r 



J Li ff (9,0',X,Z) (3»(9',Z,Y) — dZ = a vI (9,X,Y) (4.38) 



where 



L eff (9,9',X,Z) = Li ff (9, 9', X, Z) + ^(9,9', X,Z) 
a u {9,X,Y) = a uI (9,X,Y) + a" H (9,X,Y) 

and 

Ll ff (9,9',X,Z) = k5'(9-9')5(X - Z) + e8(9 - 9')5'{X - Z) - 

+ Me, M, x) L(x) + -£ <?q m (x)) H fSlf + 



Z>1 



+ 



<Sl7f(X) 



= kS'{9-9')S(X - Z) + eS(9-9')S'(X - Z) - 



VJjjW^Z) = -- J dWM8,[V],X)sgn(X -W)( R (W,Z,9')- 

a"'(0,x,y) = £e-A$(0,[u],x)«5«(x-y) 

= lf°°dWM8,[U},X)sgn(X-W) ^(W) + ^e 2 % l) (W) ) j S'(W — Y) — 



+ 
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dWij)o(6, [U], X) sgn(X -W)J2 eS K \s){ W ) s(s) ( W ~ Y ) 



We have to prove now that the solutions (3 U (9, Z, Y) satisfy in fact system (I4.36[) and 
have the form represented by fl4.32p - fl4.33p . So, let us discuss first the resolvability of 
system f)4.38p . According to relations f l4.34p -( l4~35j) we can write the main part (in e) of 
the operator L e ff in the form: 

L{ fm (6, 9', X,Z) = k 5' (9 - 9') 5(X-Z) + $^ (9, U(X)) u%(X) $ e > (&', U(X)) 5(X - Z) 

The operator L^j^ gives a set of independent operators at different X where the 
operator 

k5'(9 - 9') + ^(9,V(X))^(X)^(9',U(X)) 

has at every X exactly two linearly independent left eigen- vectors on the space of periodic 
functions in 9 

Vl (9,X) = 1 , V2 (9,X) = $(0,U(X)) 

corresponding to the zero eigen-values. 

The vectors i]i(9, X) 5(V — X) and 772(6', X) 5(V — X) give the main parts of the vectors 



6U°(V) 5U l (V) 

) \ K = 5(V-X) , ° ,} > \k = ip(9,[lJ},X)6(V-X) (4.39) 



5ip{9,Xy v ' ' 5ip(9,Xy 

which are the left eigen- vectors of the operator L^j- corresponding to the zero eigen-values. 

It's not difficult to see now that the orthogonality of the values {g(9, X),U U (Y)} to 
vectors (14.391) on K, implies the orthogonality of a uI (9, X, Y) to the same vectors. So, we 
get that system (I4.38P is a compatible system which can be resolved recursively in all the 
orders of e. The right-hand part of system (I4.38P has the form analogous to (I4.32p - fl4.33p 
so it's not difficult to see that all the (3 U (9, Z,Y) have the necessary form in this case. 
Using also the fact a uI (9, X, Y) = 0(e) we get that the solutions (3 V (9, Z, Y) have exactly 
the required form f l4.32p ~ p.33p being written as formal series in e. Besides that, condition 
(I4.37P can be also derived from a not complicated analysis of system f l4.38j) by use of the 
same left eigen-vectors of U e ^ corresponding to the zero eigen-values. 

Finally, let us prove the relation 

/ L eff (9, 9', X, Z) P(9', Z, Y)—dZ = a"(9, X, Y) (4.40) 

for the p"(6, Z,Y) found from fl438|) . 

Let us note that the difference in the images of the operators L e ff and L^j- for our 
(3 U (9,Z,Y) is proportional to ipe{9, [U],X) at every (X,Y). The same is also valid for 
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a vII (6, X, Y) which is the difference between a u (9,X,Y) and a vI (9,X,Y). It's not dif- 
ficult to check then that relation (I4.40p follows from (I4.38P and the orthogonality of the 
values 

/+oo p2w jni 
/ L% f (0, 9', X, Z) F{9', Z, Y)—dZ- a M {9, X, Y) 
-oo JO 271 

to the vectors SU 2 (V)/5(p(9, X) which takes place for our (3 U (9, Z, Y). 

Using formula (14.31 j) we can claim now that the restricted Poisson bracket 
{U V (X), U^(Y)} D has exactly form (jP3) . 

Let us recall now the the functionals I v = f_ °° U v (X)dX leave invariant the subman- 
ifold K. as was proved in Theorem 4.1. This means in particular that the flows generated 
by I v on /C coincide with their flows generated in the Dirac Poisson structure on this 
submanifold. Thus, we obtain that the functionals l v play the role of the Hamiltonian 
functions for the higher deformed Whitham systems (" 14 . 5 [) and, in particular, the func- 
tional I 2 plays the role of the Hamiltonian function for the deformed Whitham system 
(13.1 6ft after the restriction on /C. In the same way the functionals 1° and I 1 play the role 
of the annihilator and the momentum functional for the restricted Gardner - Zakharov - 
Faddeev bracket respectively. 

At last, let us say that the proof of the Theorem for the case of the Magri bracket 
repeats completely the proof for the Gardner - Zakharov - Faddeev case. 

Theorem is proved. 

Remark 4.1. 

It's not difficult to see that the main (~ e) term of the Dirac bracket {U U (X), U^{Y)}d 
on /C coincides with the Dubrovin - Novikov bracket for the Whitham system given by 
the "averaging procedure" described above. The Dubrovin - Novikov bracket obtained 
from the Gardner - Zakharov - Faddeev bracket and the Magri bracket respectively are 
compatible with each other and give a bi-Hamiltonian structure for the pure Whitham 
system for KdV. However, we can not claim here the same property for the case of the 
Dirac brackets obtained as the restrictions of the Gardner - Zakharov - Faddeev bracket 
and the Magri bracket on /C since the Dirac procedure does not preserve the compatibility 
of the brackets in general case. 

Remark 4.2. 

It's not difficult to see that the functional 

/+oo 
k(X) dX 
-oo 

plays the role of annihilator for the restricted Poisson brackets both in the cases of the 
Gardner - Zakharov - Faddeev bracket and the Magri bracket. This circumstance is 
connected with the conservation of the value S'(+oo) — S(— oo) by the flows generated by 
the "modified" functionals 
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/+oo 
U"(X) q(X) dX 
-oo 

with q(X) having compact support and has a general character for the restricted field- 
theoretical Poisson brackets. 



5 Some remarks on the averaging of the Lagrangian 
structures. 

At the end let us discuss also the averaging of Lagrangian functions for the deformed 
Whitham systems. We will restrict ourselves here only to the situation of the local La- 
grangian functions which was considered first by Whitham ( [7CTI [7T| [72]) in connection 
with the pure Whitham approach. 

As it is well-known the Gardner - Zakharov - Faddeev bracket corresponds to the local 
Lagrangian formalism of the KdV equation (13.181) having the form 
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Sv(X,T) 
(where ip = evx) which gives 



e .3 i 2 2 



-v x v T - ~v x + e v xx 



dXdT = (5.1) 



v X t + tv x vxx + e vxxxx = (5.2) 

We introduce also the Whitham pseudo-phase S(X, T) and look for the solution of 
(15. 2j) having the form 

v(9,X,T) = vw(?£^ + 0,X,T\ + ^11 (5.3) 
We require now that V^ tot \9, X, T) is a periodic function in 9 having the form 

V^ tot \9,X,T) = ^2v ik) {9,X,T) (5.4) 

fc>0 

where all the functions V^){9, X, T) are local functionals of 

(k = Sx, St, n = T, x , kx, Stx, n>x, kxx, Stxx, n xx, • • • ) 
having degree k according to the gradation rule: 

1) All the functions f(k, St,tl) have degree 0; 

2) The derivatives kkx, Sxkx, "^kx have degree k\ 

3) The degree of the product of functions having certain degrees is equal to the sum 
of their degrees. 
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According to the normalization of H(X,T) we put the conditions 



27r in 

V (k) (9,X,T) — = (5.5) 
o Zn 

for aRV (k) {9,X,T). 

Let us note that the choice of the parameters (k, St, n) instead of (k, A, n) is more 
convenient for the consideration of the Lagrangian structures in our approach. We remind 
also that the expression for St is given by relation (I3.15p . 

Easy to see that the form (15.31) gives the form of <fi(9, X, T) we consider and all the 
functions V( k )(9, X, T) are uniquely defined by the terms of series (j3.17p . Indeed, let us 
first re-expand series (I3.17P according to the new gradation rule, i.e. 



where all $^ have degree I according to the rules formulated above. 
Then we have 

kV [l)e (9,X,T) + V {l „ 1)x {9,X,T) = &^(9,X,T) , 1>1 
which defines uniquely all V^(9,X, T) view normalization rule (15. 5p and we have 



■2tt 



d9 

£x(X,T) = J <j>{9, X,T) - 
Finally, we can substitute series (15. 3p in the Lagrangian principle 

6 J J J C(9,X,T) T^dXdT 

with the Lagrangian density 

C = -S X S T ( v} tot) Y - SxSt - Isl ( V^Y - H x Sl ( v} tot) Y - Ie» + Si ( V^Y + 



xo T yv 6 j -z, x iu T --o x yv e j -n x o x yv e j -^ x ^o x yv ee j 



+ e {-S x vt t] V T tot) - S T vi tot) V { x 0t) - Y.tV^ - Y, x V T tot) - 
-S 2 X (v 9 {tot) y Vt l) - 2Z x S x V e {tot) V x tot) - Y? x V x tot) + S x V^ ot) V^ + 



_ e 2 (_y(**) V W _ Sx y(tot) (yitot)^ _ ^ fy(tot)^ + AS 2 x fy(M)^ + 
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+ s xx (yt t] ) 2 + s\v^vtx ] + ±s x s xx vt t] v^ + z xx ) + 



+ e 



3 / fiA tot )\ 3 j_ a Q \A tat )\A tat ) j_ no ir( tot )iAtot) 9V xAtoi) \ 4 / ' T/ (^o^)^ 2 
V x J +4^^ K xx + ZbxxV e V xx + ZL, XX V XX l+e 



3 

The averaged Lagrangian function 



(£)(X,T) = JC(9,X,T) 
Jo 



2;T (10 
2^ 



can be also represented in the graded form with respect to the parameters {k = Sx, St, n = 
Ex) and the Lagrangian equations 

+00 ^+00 r p+co r+00 

(C)(X,T)dXdT , ^ / / (£}(X,T)dXdT (5.6) 



SSiX^J^ x /v ' ' ' 5S(X,T)7_ 



oo ./ — oo 



give a system equivalent to (I3.16p . 

The Hamiltonian formalism in the parameters (k, St, n) can be written using La- 
grangian formalism (I5.6p . We get then the Poisson bracket in the canonical form: 

{n(X),n(Y)} = 8'{X-Y) , {*(*), J(Y)} = 8'{X — Y) (5.7) 
where J(X) is given by the graded expression 

T(x) m _ _d_<m_ d(c) v g d(c) 

{ ' dS T OX dS TX OX 2 dS T xx 4^ ^ dS Ts x 

s>0 

The Hamiltonian function is given also by the standard expression 

/+oo 
(-n(X)Z T (X) + J{X)S T {X) - (£)) dX 
-oo 

Hamiltonian structure (15. 7p is given in fact by the restriction of the Symplectic struc- 
ture corresponding to the Gardner - Zakharov - Faddeev bracket to the submanifold /C 
and so gives the canonical form of the restriction of this bracket considered in Theorem 
4.2. The functional 

+ 0O 

J(X) dX 

-oo 

gives the third annihilator of the restricted Gardner - Zakharov - Faddeev bracket, so 
we have here the complete set of the canonical variables. Finally, let us say that the 
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functionals J(X), H should be re-expanded in the graded form corresponding to the 
variables (k, A, n) using relation (13.151) to come back to our initial gradation rules. 

Let us say at the end that we believe that the averaging of the Symplectic structure is 
also possible in the case of the Magri bracket. However, the Symplectic form is much more 
complicated in this case and no procedure of such kind of Symplectic forms is known by 
now. Let us mention also, that the procedure of the restriction of Poisson and Symplectic 
structures must be also generalized to the so-called weakly non-local structures, however, 
we will not discuss here this questions. 

The author is grateful to Prof. B.A. Dubrovin for many fruitful discussions. 

The work was supported by the grant of President of Russian Federation (MD-4903.2008.2) 
and RFBR (09-01-92442-KE-a, 09-01- 12 148-ofi-m). 

References 

[1] M.J. Ablowitz, D.J. Benney., The evolution of multi-phase modes for nonlinear dis- 
persive waves, Stud. Appl. Math. 49 (1970), 225-238. 

[2] V.V. Avilov, S.P. Novikov., Evolution of the Whitham zone in KdV theory, Soviet 
Phys. Dokl. 32 (1987), 366-368. 

[3] V.V. Avilov, S.P. Novikov., Evolution of the Whitham zone in KdV theory, Soviet 
Phys. Dokl. 32 (1987), 564-566. 

[4] E.D.Belokolos, V.Z.Enolskii., Reduction of abelian functions and completely inte- 
grable equations, J. Math. Sci, vol. 106 no. 6, pp. 3395-3486, 2001 (Part I) and vol. 
108 no. 3, pp 295-374, 2002 (Part II). 

[5] F. Cooper, A. Khare, U. Sukhatme., Supersymmetry in Quantum Mechanics., World 
Scientific, Singapore 2001. 

[6] S.Yu. Dobrokhotov and V.P.Maslov., Konechnozonnye pochti periodicheskie resh- 
eniya v WKB priblizhenii., Itogi Nauki, Ser. Matem. 1980, T. 15, 3-94, (in Russian)., 
Translation: S.Yu. Dobrokhotov and V.P.Maslov., Finite-Gap Almost Periodic Solu- 
tions in the WKB Approximation. J. Soviet. Math., 1980, V. 15, 1433-1487. 

[7] S. Yu. Dobrokhotov., "Resonances in asymptotic solutions of the Cauchy problem for 
the Schrodinger equation with rapidly oscillating finite-zone potential" , Mathematical 
Notes, 44:3 (1988), 656-668. 

9 Lct us note also that relation (|3.15[) can not be denned from Lagrangian function and should be 
defined separately from the asymptotic procedure. 



49 



S. Yu. Dobrokhotov. , "Resonance correction to the adiabatically perturbed finite- 
zone almost periodic solution of the Korteweg - de Vries equation", Mat. Zametki, 
44:4 (1988), 551-555. 

B.A. Dubrovin., Functionals of the Peierls - Frohlich type and the variational principle 
for the Whitham equations, . Amer. Math. Soc. Transl. (2) 179 (1997), 35-44. 

B. Dubrovin, On Hamiltonian perturbations of hyperbolic systems of conservation 
laws, II: universality of critical behaviour, Comm. Math. Phys. 267 (2006) 117 - 139. 

B.A. Dubrovin and S.P.Novikov., Hamiltonian formalism of one-dimensional systems 
of hydrodynamic type and the Bogolyubov - Whitham averaging method, Soviet 
Math. Dokl, Vol. 27, (1983) No. 3, 665-669. 

B.A. Dubrovin and S.P. Novikov., Hydrodynamics of weakly deformed soliton lat- 
tices. Differential geometry and Hamiltonian theory, Russian Math. Survey, 44 : 6 
(1989), 35-124. 

B. A. Dubrovin and S.P.Novikov., Hydrodynamics of soliton lattices, Sov. Sci. Rev. 

C, Math. Phys., 1993, V.9. part 4. P. 1-136. 

B.A. Dubrovin., "Integrable systems in topological field theory", Nucl. Phys., B379 
(1992), 627-689. 

B.A. Dubrovin., Integrable Systems and Classification of 2-dimensional Topological 
Field Theories, ArXiv: |hep-th/9209"040| 

B.A. Dubrovin., Geometry of 2d topological field theories, ArXiv: hep-th/9407018 



B.A. Dubrovin., "Flat pencils of metrics and Frobenius manifolds" 



ArXiv: math. DG/9803 106, In: Proceedings of 1997 Taniguchi Symposium 



"Integrable Systems and Algebraic Geometry", editors M.-H.Saito, Y.Shimizu and 
K.Ueno, 47-72. World Scientific, 1998. 

B.A. Dubrovin., "Geometry and analytic theory of Frobenius manifolds", ArXiv: 
math. AG/98070341 



B.A. Dubrovin, Y.Zhang., Bihamiltonian Hierarchies in 2D Topological Field Theory 
At One-Loop Approximation, Commun. Math. Phys. 198 (1998), 311-361. 

B.A. Dubrovin, Y.Zhang., Normal forms of hierarchies of integrable PDEs, Frobenius 
manifolds and Gromov-Witten invariants., ArXiv: math.DG/0108160 

B.A. Dubrovin, Y.Zhang., Virasoro Symmetries of the Extended Toda Hierarchy, 
ArXiv: |math.DG/0308152| 

Boris Dubrovin, Si-Qi Liu, Youjin Zhang., "On Hamiltonian perturbations of hyper- 
bolic systems of conservation laws", ArXiv: math.DG/0410027 



50 



Boris Dubrovin, Youjin Zhang, Dafeng Zuo., "Extended affine Weyl groups and 
Frobenius manifolds - II", ArXiv: math.DG/0502365 

G.A. El, A.L. Krylov, S. Venakides., Unified approach to KdV modulations, Comm. 
Pure Appl. Math, Comm. Pure Appl. Math, 54 : 10 (2001), 1243-1270. 

E.V. Ferapontov., Differential geometry of nonlocal Hamiltonian operators of hydro- 
dynamic type, Functional Anal, and Its Applications, Vol. 25, No. 3 (1991), 195-204. 

E.V. Ferapontov., Dirac reduction of the Hamiltonian operator S 1 ^-^ to a subman- 
ifold of the Euclidean space with flat normal connection, Functional Anal, and Its 
Applications, Vol. 26, No. 4 (1992), 298-300. 

E.V. Ferapontov., Nonlocal matrix Hamiltonian operators. Differential geometry and 
applications, Theor. and Math. Phys., Vol. 91, No. 3 (1992), 642-649. 

E.V. Ferapontov., Nonlocal Hamiltonian operators of hydrodynamic type: differential 
geometry and applications, Amer. Math. Soc. Transi, (2), 170 (1995), 33-58. 

Flaschka H., Forest M.G., McLaughlin D.W., Multiphase averaging and the inverse 
spectral solution of the Korteweg - de Vries equation, Comm. Pure Appl. Math., - 
1980.- Vol. 33, no. 6, 739-784. 

T. Grava., Self-similar asymptotic solutions of the Whitham equations, Russian 
Math. Survey, 54 : 2 (1999), 169-170. 

T. Grava., Existence of a global solution of the Whitham equations, Theor. Math. 
Physics 122 : 1 (2000), 46-57. 

T. Grava., From the solution of the Tsarev system to the solution of the Whitham 
equations, Mathematical Physics, Analysis and Geometry 4 : 1, (2001), 65-96. 

T. Grava., Riemann-Hilbert problem for the small dispersion limit of the KdV equa- 
tion and linear overdetermined systems of Euler-Poisson-Darboux type, Comm. Pure 
Appl. Math, 55 : 4 (2002), 395-430. 

T. Grava, F.-R. Tian., The generation, propogation and extinction of multiphases in 
the KdV zero dispersion limit, Comm. Pure Appl. Math, 55 : 12 (2002), 1569-1639. 

N.I.Grinberg., On Poisson brackets of hydrodynamic type with a degenerate metric, 
Russian Math. Surveys, 40:4 (1985), 231-232. 

A.V. Gurevich, A.L. Krylov, and G.A. El., Riemann wave breaking in dispersive 
hydrodynamics, J FTP Letters 54 (1991), 102-107. 

A.V. Gurevich, A.L. Krylov, and G.A. EL, Evolution of a Riemann wave in dispersive 
hydrodynamics, Sov. Phys. JETP 74 (1992), 957-962. 



51 



A.V. Gurevich, L.P. Pitaevskii., Decay of initial discontinuity in the Korteweg - de 
Vries equation, JETP Letters 17 (1973), 193-195. 

A.V. Gurevich, L.P. Pitaevskii., Nonstationary structure of a collisionless shock 
waves, Sov. Phys. JETP 38 (1974), 291-297. 

A.V. Gurevich, L.P. Pitaevskii., Averaged description of waves in the Korteweg - de 
Vries - Burgers equation, Soviet Phys. JETP 66 (1987), 490-495. 

W.D. Hayes., Group velocity and non-linear dispersive wave propagation, Proc. Royal 
Soc. London Ser. A 332 (1973), 199-221. 

I.M. Krichever., The averaging method for two-dimensional integrable equations, 
Functional Anal. Appl. 22 (1988), 200-213. 

I.M. Krichever., Perturbation theory in periodic problems for two-dimensional inte- 
grable systems, Sov. Sci. Rev. Section C9 (1992). 

I.M. Krichever, S.P. Novikov., Evolution of the Whitham zone in the Korteweg - de 
Vries theory, Soviet Phys. Dokl. 32 (1987), 564-566. 

P.D. Lax, CD. Levermore., The small dispersion limit for the Korteweg - de Vries 
equation I, II, and III. Comm. Pure Appl. Math., 36 (1983), 253-290, 571-593, 809- 
830. 

P.D. Lax, CD. Levermore, S. Venakides., The generation and propagation of oscil- 
lations in dispersive IVPs and their limiting behavior, Important developments in 
soliton theory 1980-1990, 205-241, Springer Series in Nonlinear Dynamics. Springer, 
Berlin (1993). 

Si-Qi Liu, Youjin Zhang., Deformations of Semisimple Bihamiltonian Structures of 
Hydrodynamic Type, ArXiv: |math.DG/0 405146 

Si-Qi Liu, Youjin Zhang., On the Quasitriviality of Deformations of Bihamiltonian 
Structures of Hydrodynamic Type, ArXiv: |math.DG/ 0406626 

P. Lorenzoni., Deformations of bihamiltonian structures of hydrodynamic type, J. 
Geom. Phys. 44 (2002), 331-371. 

Luke J.C, A perturbation method for nonlinear dispersive wave problems, Proc. Roy. 
Soc. London Ser. A, 292, No. 1430, 403-412 (1966). 

A.Ya.Maltsev., The conservation of the Hamiltonian structures in Whitham's method 
of averaging, Izvestiya, Mathematics 63:6 (1999), 1171-1201. 

A.Ya.Maltsev., The averaging of non-local Hamiltonian structures in Whitham's 
method. // solv-int/9910011, International Journal of Mathematics and Mathemat- 
ical Sciences, 30:7 (2002) 399-434. 



52 



A.Ya.Maltsev, S.P. Novikov. -"On the local systems hamiltonian in the weakly non- 
local Poisson brackets." ArXiv: |nlin.SI/0006030~[ Physica D 156 (2001) 53-80. 

A.Ya.Maltsev. "Weakly- nonlocal Symplectic Structures, Whitham method, 
and weakly-nonlocal Symplectic Structures of Hydrodynamic Type." Arxiv: 
|nlin.SI/0405060 , J. Phys. A: Math. Gen. 38 (3) 637-682. 

A.Ya. Maltsev., Whitham systems and deformations, Journ. Math. Phys. 47, (2006). 

A. Ya. Maltsev. The deformation of the Whitham systems in the almost linear case, 
arXiv: 0709.4618 

O.I. Mokhov and E.V. Ferapontov., Nonlocal Hamiltonian operators of hydrodynamic 
type associated with constant curvature metrics, Russian Math. Surveys, 45:3 (1990), 
218-219. 

A. C. Newell. Solitons in mathematics and physics. Society for Industrial and Applied 
Mathematics (1985). 

S.P. Novikov., The geometry of conservative systems of hydrodynamic type. The 
method of averaging for field-theoretical systems, Russian Math. Surveys. 40 : 4 
(1985), 85-98. 

M.V.Pavlov., Elliptic coordinates and multi-Hamiltonian structures of systems of 
hydrodynamic type., Russian Acad. Sci. Dokl. Math. Vol. 59 (1995), No. 3, 374-377. 

M.V.Pavlov., Multi-Hamiltonian structures of the Whitham equations, Russian 
Acad. Sci. Doklady Math., Vol. 50 (1995) No.2, 220-223. 

G.V. Potemin, Algebraic-geometric construction of selfsimilar solutions of Whitham 
equations, Russian Math. Surveys 43 : 5 (1988), 252-253. 

F.-R. Tian., Oscillations of the zero dispersion limit of the Korteweg - de Vries 
equation, Comm. Pure Appl. Math 46 (1993), 1093-1129. 

F.-R. Tian., The initial value problem for the Whitham averaged system, Comm. 
Math. Phys 166 (1994), 79-115. 

S.P.Tsarev., On Poisson brackets and one-dimensional Hamiltonian systems of hy- 
drodynamic type, Soviet Math. Dokl, Vol. 31 (1985), No. 3, 488-491. 

S. Venakides., The zero dispersion limit of the KdV equation with non-trivial reflec- 
tion coefficient, Comm. Pure Appl. Math., 38 (1985), 125-155. 

S. Venakides., The generation of modulated wavetrains in the solution of the Ko- 
rteweg - de Vries equation, Comm. Pure Appl. Math., 38 (1985), 883-909. 



53 



[68] S. Venakides., The zero dispersion limit of the periodic KdV equation, Trans. Amer. 
Math. Soc. 301 (1987), 189-226. 

[69] S. Venakides., The Korteweg de Vries equations with small dispersion: higher order 
Las - Levermore theory, Comm. Pure Appl. Math. 43 (1990), 335-361. 

[70] G. Whitham, A general approach to linear and non-linear dispersive waves using a 
Lagrangian, J. Fluid Mech. 22 (1965), 273-283. 

[71] G. Whitham, Non-linear dispersive waves, Proc. Royal Soc. London Ser. A 139 
(1965), 283-291. 

[72] G. Whitham, Linear and Nonlinear Waves. Wiley, New York (1974). 

[73] V.E. Zakharov, S.V. Manakov, S.P. Novikov, L.P. Pitaevskii., Teoriya solitonov. 
Metod obratnoi zadachi. Nauka, Moscow 1980. (ed. S. P. Novikov) (in Russian)., 
Translation: S.P. Novikov, S.V. Manakov, L.P. Pitaevskii, and V.E. Zakharov., The- 
ory of solitons. The inverse scattering method., Plemun, New York 1984. 



54 



